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Foreword 


Proper  support  of  Mine  Warfare,  Mine  Countermeasures,  Amphibioas 
Warfare,  Special  Warfare,  and  Antisubmarine  Warfare  requires 
techniques  and  equipment  to  characterize  the  surface  and  near-surface 
sea-floor  sediments  in  real-time  from  vessels  while  underway. 

The  Naval  Ocean  Research  and  Development  Activity  (NORDA)  is 
investigating  and  developing  acoustic,  electric  resistivity,  (DC)  and  elec¬ 
tromagnetic  (EM)  technologies  for  application  to  the  real-time,  under¬ 
way  sediment  classification  task.  This  report  shows  that  DC  and  EM 
technologies  have  potential  practical  application  in  the  underway  sedi¬ 
ment  classification  task;  work  continues  in  clarifying  the  limits  of  applica¬ 
tion  of  this  ':^chnology  in  the  Navy  mission  execution. 


A.  C.  EsaUy  Captain,  USN 
Commanding  Officer,  NORDA 


Executive  summary 


The  initial  results  of  a  technical  feasibility  study  for  using  direct  current 
(DC)  electric  resistivity,  electromagnetic  (EM) ,  and  induced  polarization 
(IP)  methods  to  classify  surficial  sea-floor  sediments  are  reported, -herein. 
Analytical  work  performed  in  this  study  and  a  review  of  related  work 
show  that  both  the  DC  and  EM  resistivity  techniques  can  provide 
electrical  resistivity  measurements  of  sufficient  resolution  to  classify 
sediments  according  to  type  (mud.  clay,  silt,  sand,  gravel,  etc.).  Both  DC 
and  EM  antennas  must  ^  close  to  the  sea-floor  surface  (within  10  m)  so 
that  surficial  sediment  resistivities  can  be  sufficiently  resolved  to  permit 
sediment  classification.  To  resolve  the  resistivity  of  a  1-m-thick  surficial 
sea-floor  sediment  layer,  the  DC  antennas  must  be  deployed  less  than  1 
m  above  the  sea  floor.  Both  DC  and  EM  methods  can  adequately  resolve 
the  resistivity  of  a  6-m-thick  surficial  layer  at  a  sundoff  distance  of  10 
m.  Noise  tests  with  an  inverted  Schlumberger  array  towed  in  the  water 
column  and  on  the  sea  floor  show  noise  levels  to  be  within  acceptable 
limits.  The  level  of  development  of  the  IP  method  for  clay  mineral 
detection  is  reported;  i  Js  method  shows  considerable  promise  as  a 
contributor  to  a  remote  underway  sediment  classification  system. 
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Rapid  underway  sediment  classiflcation 
by  electrical  methods  _ 


Introduction 

Purpose 

This  report  f>resents  the  results  of  work  conducted 
during  fiscal  year  1986  in  the  evaluation  of  direct 
current  (DC)  and  electromagnetic  (EM)  technologies 
for  use  in  the  geoacoustic/geotechnical  characterization 
of  sea-floor  sediments. 

Background 

The  Navy  requires  more  rapid  techniques  and 
equipment  to  characterize  the  surface  and  near-surface 
sea-'floor  sediments  to  properly  support  Antisubmarine 
Warfare  (ASW),  Mine  Warfare  (MIW),  Mine  Counter¬ 
measures  (MCM),  and  Amphibious  Forces.  Present 
techniques  for  this  sediment  charaaerization  include 
measurement  of  properties  on  sediment  cores  eitl  er  on 
a  ship  or  in  the  laboratory,  in  situ  measurement  of 
sediment  properties  using  instrumented  free-fall 
penetrometera  and  instrumented  probes  from  bottom¬ 
resting  platforms,  and  calculation  of  properties  through 
inversion  of  geophysical  reflection  and  refraction  data. 

The  first  two  methods  are  labor  intensive  and/or 
expensive  and,  more  importantly,  both  corers  and  in 
situ  penetrometers/probes  gather  data  only  at  the 
specific  location  sampled,  allowing  deviations  between 
sample  locations  to  pass  undetected.  The  third  method 
samples  nearly  continuous  along  a  trackline,  but  the 
data  measured  are  generally  usable  only  to  the 
Antisubmarine  Warfare  community;  even  then,  resolu¬ 
tion  in  the  surficial  sediments  is  less  than  desired. 
Techniques/equipment  that  can  be  applied  from  an 
underway  vessel  are  highly  desirable  or  mandatory  to 
meet  Navy  requirements.  The  need  for  an  underway 
system  leads  to  the  selection  of  remote  sensing  systems, 
using  acoustic,  DC  resistivity  and/or  EM  induction 
sensing  technologies,  with  the  aid  of  possible  ground- 
truthir  using  penetrometers  or  sediment  cores. 

An  acoustic  sediment  classification  system  has  been 
developed  and  is  being  improved  anr’i  adapted  to 
specialized  applications  by  the  Naval  Ocesn  Research 
and  Development  Activity  (NORDA)  for  various 


sponsors.  An  evaluation  of  the  accuracy  of  sediment 
type  and  properties  predicted  by  the  acoustic  system 
over  sea-fioor  sediments  ranging  from  soft  muds  to 
dense  sands  (Lambert,  1987)  indicates  the  acoustic 
system  does  properly  classify  surficial  continental  shelf 
sedimenu  about  90%  of  the  time,  when  ground  truth 
dau  are  available  to  constrain  the  sediment  type  and 
property  predictions.  With  the  addition  of  a  second, 
remote  sensing  technology,  i.e.,  one  or  both  of  the 
electrical  technologies,  it  appears  possible  to  improve 
the  accuracy  of  the  remote  classifications  to  provide 
correct  predictiorts  95-98%  of  the  time,  with  a 
reduction  in  the  number  oi  ground  truth  data  required. 
In  addition,  the  electrical  techrtiques  are  capable  of 
performing  sediment  classificatiorts  where  the  acoustic 
methods  cannot  perform  adequately:  (1)  in  classifying  a 
sediment  of  low  impedance  (soft  mud)  beneath  a  layer 
of  higher  acoustic  impedance  (sand),  (2)  in  properly 
classifying  gassy  sediments,  (3)  in  classifying  sediments 
below  the  subseafloor  design  classification  depth  of  the 
present  acoustic  system  of  6  m  in  soft  clays  and  1  m  in 
sands,  and  (4)  in  measuring  the  sediment  thickness  over 
rock  of  low  permeability. 

Expendable  geotechnical  penetromet.  i  systems  are 
available  for  obtaining  ground-truth  data  from  a  vessel 
underway  (Beard,  1984),  and  a  more-economical, 
shallow-water,  recoverable  penetrometer  system  can  be 
assembled  from  existing  hardware  and  technology 
(Valent,  1986).  Thus  the  ground-truthing  portion  of  the 
remote  underway  system  is  available. 

Thus,  the  electrical  techniques  are  the  only  elements 
of  the  rapid  sea-floor  sediment  classification  system  not 
near-operational  at  this  time.  Direct  current  resistivity 
and  electromagnetic  techniques  are  highly  developed  for 
terrestrial  use  for  locating  and  delineating  bodies  of 
certain  mineral  deposits  and  deposits  of  water  bearing 
gravels  and  sands.  In  the  oceans.  DC  resistivity 
techniques  have  been  used  on  the  continental  shelf  in 
the  Arctic  to  measure  the  unfrozen  sediment  thickness 
over  permafrost  and  to  measure  the  thickness  of  the 
permafrost  (Corwin,  1983).  Electromagnetic  tech¬ 
niques  employed  from  aircraft-towed  sensors  in  a  bird 
have  been  used  for  the  measurement  of  sea  ice  thick- 
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ness  and  water  depth  in  the  near  shore,  and  may  prove 
capable  of  classifying  sediments  in  shallow  ocean  depths 
(i.e.,  in  the  near  shore)  (Won  and  Smits,  198S).  The 
concept  of  mounting  the  EM  source  and  receiver  in  a 
towfish  permits  the  classifier  hardware  to  be  flown 
closer  to  the  sea-floor  sediment  being  scrutinized,  and 
thereby  improves  the  reliability  of  the  EM  system 
sediment  classification  and  its  resolution  of  sediment 
layers.  The  viability  of  such  a  system  was  first  proposed 
by  Hulbert  et  al.  (1982). 

Approach 

The  technical  feasibility  of  using  DC  and  EM 
techniques  for  classifying  sea-floor  material  has  been 
initiated  through  the  use  of  analytical  techniques.  The 
FY86  effort  was  a  quick-look  evaluation  of  the 
influence  of  stand-off  disunce  between  tlie  array  and 
the  sea  floor:  the  capability  to  discriminate  between 
sediment  types;  and  the  probable  geometries,  sizes,  and 
power  requirements  of  the  arrays.  Work  on  the  DC 
technology  area  was  subcontraaed  to  Harding-Lawson 
Associates.  Novato.  California.  (Dr.  Robert  F.  Corwin). 
Dr.  Corwin  delivered  a  collection  of  measured  in  situ 
sediment  electrical  resistivities,  the  results  of  a  computer 
study  of  the  performance  of  a  DC  array  used  in  the 
sediment  classiflcation  mode,  and  the  results  obtained 
from  field  measurements  of  signal-to-noise  ratio  in  a 
full-scale  DC  array  towed  in  San  Francisco  Bay 
(Corwin.  1986).  In-house  studies  have  first  developed 
and  validated  a  computer  algorithm  to  predict  the 
performance  of  an  EM  system  with  both  source  and 
receiver  towed  in  the  seawater  column.  This  algorithm 
has  been  exercised  to  produce  some  preliminary  results 
of  EM  system  performance  reported  herein  This  report 
presents  the  results  of  these  FY86  efforts,  suiting  with  a 
brief  discussion  of  electrical  resistivity  in  sediments, 
followed  by  a  description  of  the  DC  and  EM  techniques, 
and  finally  the  performance  results  provided  by  the 
computer  analyses. 

Electrical  resistivity  of 
sea-floor  materials 

The  phenomena 

Electrical  conduction  in  most  indurated  sediments 
and  rocks  and  in  most  cohesionless  sedimenu  is 
essentially  electrolytic  via  the  material  pore  fluid, 
because  most  mineral  grains  are  insulators.  The 
seawater  pore  fluid  of  marine  rocks  and  sediments 
provides  the  ideal  medium  for  the  conduction  of  an 
electrical  current,  with  the  resistivity  of  the  formation 
inversely  related  to  the  effective  porosity  and  the  degree 
of  saturation.  The  effective  porosity  is  a  fonction  of  the 
measured  porosity,  i.e-,  the  ratio  of  the  volume  of  voids 
to  the  toul  volume,  and  of  the  grain  size  and  shape,  and 


pore  size  and  shape.  When  sediments  and  rocks  con^-ain 
a  measurable  portion  of  clay  minerals  and/or  metallic 
sulfides  and  magnetite,  the  sediment  matrix  materials  do 
Uke  part  in  the  electrical  conduction  process.  With  clay 
minerals,  conduction  ukes  place  in  the  cation  double 
layer  surrounding  the  actual  miners!  particle. 

The  resistivities  of  terrestrial  geologic  materials  cover 
a  very  wide  range,  from  less  than  IDm  to  greater  than 
10*fiin  (Griffiths  and  King,  1965;  Keller  and 
Frischknecht,  1966).  The  resistivities  of  sea-floor 
sedimenu  are  much  lower,  ranging  from  O.Sflm  to 
26flm  (dau  referenced  by  Corwin,  1986),  primarily 
because  of  the  highly  conductive  seawater  pore  fluid. 

Formation  Factor 

Because  the  resistivity  of  a  geological  formation  is 
largely  dependent  on  the  resistivity  of  the  pore  fluid,  the 
discussion  of  formation  resistivities  is  simplified  by  use  of 
the  formation  factor  (F),  defined  as  (Archie,  1942); 

F  «  iy7p.».  (1) 

where  py  »  resistivity  of  the  bulk  formation,  and  p^,  a 
resistivity  of  the  pore  water.  For  terrestrial, 
consolidated,  relatively  low-porosity  rock,  F  is  related 
to  porosity  (n)  by  the  empirical  formula  knovm  as 
Archie’s  law: 

F  ■  an*'",  (2) 

whei'e  o,  m  are  constants  with  values  close  to  1  and  2. 
resptictively  (Keller  and  Frischknecht,  1966). 

I'he  resulu  of  thousands  of  measurements  verify  the 
direct  relationship  between  n  and  F  for  terrestrial  rocks 
of  low  clay  content.  However,  because  clay  minerals 
readily  conduct  an  electric  current  in  their  surface 
dcsible  layer,  the  measured  electrical  resistances  (py) 
and  formation  factors  (F)  for  clay-rich  rocks  and 
sedimenu  may  be  significantly  lower  than  predicted  by 
Archie's  law  for  a  given  porosity. 

Typical  values  of  formation  factor 

Because  the  salinity  of  seawater  is  greater  than  that  of 
most  terrestrial  sediment  pore  water,  the  absolute 
resistivity  and  formation  factors  jf  marine  sediments 
and  rocks  are  generally  considerably  lower  than  that  of 
the  same  materials  on  land.  Measured  values  of  F  for 
sea-floor  sedimenu  range  from  less  than  1  for 
near-mudline  arctic  sediments  with  elevated  salt 
concentration  in  the  pore  water  to  about  IS  for  some 
coat!,'!  sands  and  gravels  (Table  1).  Formation  factors 
greater  than  7  or  8  generally  indicate  rock  or  bonded 
permafrost  rather  than  unfrozen  sedimentary  material. 

There  is  no  unique  relationship,  applicable  world¬ 
wide,  between  resistivity  and  sea-floor  sediment 
propr^rties.  On  a  world-wide  basis,  measured  resistivities 
for  clay  may  overlap  those  for  silu,  silts  may  overlap 
sands,  and  sands  may  overlap  gravels.  However,  within 
a  survey  area,  the  relationship  between  electrical 
resistivity  and  sediment  type  is  usually  much  better 


.  Summary  of  selected  sea-floor  resistivity  studies  (from  Corwin, 
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constrcined,  and  ledimem  type  can  be  predicted  from 
the  meaturtd  resistivity. 

Generally,  a  difference  of  one  unit  in  formation  factor 
indicates  a  significant  difference  in  sediment  properties. 
Thus,  for  example,  in  a  given  survey  area  where  a  sf*ty 
sediment  has  a  formation  factor  of  2,  a  change  to  3  in 
the  measured  formation  faaor  would  Lidicate  signih- 
cantly  larger  grain  size,  decrersed  porosity,  or  increased 
degree  of  consolidation. 

Direct  current  resistivitv 
methods 

The  technique 

In  the  DC  method  for  measuring  the  electrical 
resistivity  of  the  sea-floor  materials,  a  current  is  injeaed 
into  the  sediments  or  into  the  seawater  by  means  of  two 
electrodes,  and  the  potential  difference  between  a 
second  pair  of  electrodes  placed  in  line  between  the  Hrst 
pair  is  measured  (Fig.  la  and  b).  From  the  magnitude 
of  the  ratio  of  potential  difference  to  appUed  current 
and  a  knowledge  of  the  electrode  separation,  a  quantity 
known  as  the  apparent  resistivity  can  be  calculated.  If 
the  sea-floor  material  is  homogeneous,  then  this 
measured  resistivity  is  the  same  as  the  true  ground 
resistivity,  but  in  general  it  is  a  vreighted  average  of  the 
resistivities  of  the  various  strata  through  which  the 
current  passes. 

For  marine  application,  the  inverted  Schlumberger 
array  (Fig.  Ic)  provides  the  best  combination  of  layer 
resolution,  signal-to-noise  ratio,  and  convenience  of 
operation.  In  terrestrial  applications,  the  current 
(receiver)  electrodes  are  insulled  at  the  center  of  the 
array,  the  potential  electrode  qsacing  is  expanded,  and 


Figure  I.  Example  arrangements  oj  electrodes  for  DC 
method. 


measurements  are  made  at  each  new  position  to  develop 
the  curve  of  apparent  resistivity  versus  electrode 
snpMiration.  In  marine  applications,  it  is  more  convenient 
to  deploy  the  electrodes  in  a  single  string  or  eel  behind  a 
towing  vessel,  with  the  string  towed  in  the  water  column 
or  diagged  slowly  alorg  the  sea  floor.  The  marine  system 
US3S  an  inverted  Scitiumberger  configuration  with  the 
current-transmitting  electrodes  A  and  B  deployed  at  the 
center  of  the  array  (Pig.  2a)  and  the  voitage-receiving 
electrodes  Mi.  Ms,  .  .  .  Mn  and  Ni,  Nj,  .  .  .  Nn 
placed  at  logarithmically  increasing  spacings  from  the 
center  of  the  array.  This  multiplicity  of  receiving 
eleinrodes  allows  all  the  voluge  readings  to  be  obtained 
simultaneously  from  a  single  current  pulse,  all  at 
virtually  the  ume  insunt  in  time,  so  that  quality  data 
can  be  gathered  from  a  vessel  underway. 

Optimum  array  geometries 

Values  of  the  current  and  potential  electrode 
spacings.  AB  and  MN,  are  selected  to  provide  both 
adequate  resolution  of  relatively  thin  layers  close  to  the 
sea  floor  and  adequate  depth  of  investigation  below  the 
sea  floor.  The  maximum  subseafloor  depth  of 
investigation  of  a  DC  resistivity  measurement  is  a 
complex  function  of  the  resistivity  structure  and  the 
separation  of  the  outer  electrodes  (Roy  and  Elliot, 
1961).  In  general,  with  the  electrodes  at  the  sea-floor 
surface,  the  effective  depth  of  investigation  is  about 
one-sixth  to  one-fourth  of  the  maximum  potential 
electrode  separation  (MN).  Thus,  for  the  array  depicted 
in  Figure  2a.  with  a  maximum  MN  of  63.2  m  and  array 
deployed  on  the  sea  floor,  apparent  resistivities  and 
structure  will  be  resolved  to  a  subseafloor  deptl,  of  about 
10  to  15  m. 

The  minimum  potential  electrode  spacing  is  selected 
to  resolve  the  minimum  desired  thicknesses  close  to  the 
sea  floor.  Again,  for  the  array  of  Figure  2a,  the 
minimum  potential  electrode  separation  (AB)  of  2.0  m 
will  allow  resolution  of  the  parameters  of  a  layer  of 
about  0.3  m  thickness  at  the  sea  floor. 

Accurate  interpretation  of  data  from  a  Schlumberger 
array  requires  that  the  current  electrode  spacing  (AB) 
be  very  small  as  compared  to  the  potential  elecuode 
spacing  (MN).  A  practical  minimum  spacing  is  about 
one-fourth  of  the  potential  electrode  spacing,  or 

(AB)  recommended  ®  (IM)  (MN).  (3) 

While  optimally  Jte  current  electrode  spacing  (AB) 
must  be  small  to  permit  resolution  of  thin  layers,  the 
spacing  must  be  large  to  develop  an  adequate  received 
voluge  signal  at  the  outermost  potential  electrodes.  For 
a  given  transmitted  current  l^/el  (I)  at  the  current 
electrodes,  larger  current  electrode  spacing  (AB)  yields 
larger  potential  difference  (volume)  at  the  receiving 
electrodes.  The  received  volu^;?  or  signal  at  the 
potential  electrodes  varies  as  the  inverse  of  the  square  of 
the  disunce  between  the  electrodes  (MN)^,  so  that  the 


Figure  2.  (a)  Schematic  of  towed  marine  Schiumberger  array  for  direct  current  resistivity  profiling,  and 
(b)  a  typical  plot  of  apparent  resistivity  data  versus  receiver  electrode  spacing. 


•ignals  at  large  voltage  electiode  ipacings  are  much 
tmaller  than  those  at  the  inner  electives.  If.  to  obtain 
the  required  penetration,  the  outer  potential  electrode 
spacing  (MN)  must  be  made  very  large  and  the  received 
signal  is  then  unacceptably  small,  then  a  second  set  of 
current  electrodes  at  a  larger  spacing  (AB)  can  be  used. 
Appit^riate  sequential  switching  circuitry  then  would  be 
used  to  ensure  that  readings  are  uken  from  the 
appropriate  potential  elecuode  pairs  (MN)  when 
transmitting  from  various  current  electrode  pairs  (AB). 

Thus,  selection  of  the  current  electrode  spacing  (AB) 
and  the  minimum  and  maximum  potential  electrode 
spacings  (MN)  must  be  balanced  among  requiremenu 
for  maximum  depth  of  investigation,  desired  resolution, 
available  transmitter  current  levels,  and  acceptable 
signal'to-noise  ratios  at  the  larger  current  electrode 
spacings  (MN). 

Equipment 

The  major  components  of  a  marine  DC  resistivity 
measurement  system  include  a  current  transmitter, 
transmitting  and  receiving  electrodes,  cables  connecting 
the  electrodes  to  the  transmitter  and  receiver,  and 
instrumenution  to  measure  and  record  the  transmitted 
current  »nd  received  voluges. 

Current  transmitter 

Because  the  DC  voltage  across  the  receiving 
electrodes  exhibits  a  certain  amount  of  drift  (discussed 
below),  current  is  transmitted  in  the  form  of  a  square 
wave  with  a  frequency  of  about  0.1  to  1  Hz.  By 


measuring  the  peak-to-peak  volt*  ge  across  the  receiving 
electrodes  and  dividing  by  turn,  the  effect  of  any  DC 
offset  is  automatically  removed  and  the  true  signal 
vohage  aiiq>litude  can  be  determined. 

The  current  transmitter  system  consists  of  a  power 
rupply,  a  switching  system,  and  equipment  to  measure 
and  record  the  transmitted  current.  To  provide  usable 
signals  for  a  practical  array,  the  transmitted  current  must 
be  of  the  order  of  several  amperes  (A)  or  more  and 
have  very  low  noise  level.  Transmission  of  currents  of 
this  magnitude  into  seawater  requires  a  power  supply  of 
several  hundred  vohs.  A  3.S-kW,  gasoline-powered, 
motor-generator  and  a  control  system  that  steps  up, 
rectifies,  and  filters  the  AC  output  of  the  generator, 
providing  a  continuously  variable  output  of  0  to  800 
VDC  at  up  to  10  A  has  been  used  for  this  task.  Larger 
systems  using  the  sa-'ie  concept  could  easily  be 
constructed. 

Switching  frequencies  of  about  0.1  to  1  Hz  allows 
transient  switching  effects  to  die  out  before  the  voltage 
measurements  are  made.  Higher  frequencies  may  allow 
more  lapid  dau  acquisition,  but  transient  and  inductive 
and  capacitive  effecu  may  impair  data  accuracy. 

Switching  the  current  elecuode  polarity  can  be 
accomplished  manually  or  electronically.  Manual 
switching  offers  high  reliability  and  continuous  operator 
conuol  over  the  duty  cycle,  but  is  not  practical  for 
extended  survey  lines  and  for  operating  frequencies 
greater  than  about  1  Hz.  Electronic  switching  frees  the 
operator  for  other  tasks,  but  requires  sophisticated 
high-power  circuitry. 


IVansmitting  electrodes 

Current  trtnsmitting  (AB)  electrodes  for  marine  DC 
measurements  can  be  fabricated  as  wire  wrapped 
around  the  outside  of  the  electrode  streamer  cable 
(discussed  below).  A  length  of  about  10  cm  wrapped 
around  a  cable  of  1-  to  2-cm  diameter  gives  an  effective 
contact  resisunce  of  a  few  tens  of  ohms.  Thus  a  source 
voltage  of  about  250  VDC  will  dnve  a  current  of  about  S 
to  10  A  through  such  an  electrode.  We  have  found  both 
copper  and  Monel  wire  to  be  satisfactory  electrode 
materials,  with  minimal  metal  less  even  after  many  days 
of  operation  at  10- A  levels. 

Receiving  electrodes 

The  voltage-receiving  (MN)  electrodes  must  permit 
measuring  voltage  differences  of  a  few  tenths  of  a 
millivolt  (mV)  with  minimal  intrinsic  noise  and  drift. 
Although  bare  metallic  electrodes  of  noble  metals,  such 
as  gold  or  platinum,  are  used  for  such  measurements  in 
the  laboratory,  their  noise  and  drift  levels  increase  to 
unacceptable  levels  during  field  use.  Nonpolarizing 
silver-silver  chloride  (Ag-AgCl)  electrodes  perform  well 
for  marine  measurements  (Corwin  and  Conti,  1973) 
These  electrodes  consist  of  a  chloridized  silver  wire 
immersed  in  a  conductive  gel  contained  in  a  plastic 
housing.  Contact  between  the  seawater  and  the  gel  is 
through  a  porous  ceramic  junction.  Thus  the  wire 
electrode  element  is  isolated  from  direct  contact  with 
the  seawater,  with  a  corresponding  reduction  in  noise 
generated  by  water  turbulence  and  water  temperature 
and  salinity  variations.  Intrinsic  noise  levels  for  such 
electrodes  deployed  stationary  on  the  sea  floor  are  of 
the  order  of  a  few  microvolts  (uV),  and  long-term  drift 
is  less  than  a  few  millivolts  per  day. 

Electrodes  of  this  type  are  available  commercially 
(they  are  used  as  reference  elecuodes  for  pH  and  other 
electrochemical  measurements),  or  can  be  fabricated  to 
meet  special  requirements. 

Connecting  cables 

For  a  marine  DC  resistivity  system  all  the  transmitting 
and  receiving  electrodes  are  fabricated  into  a  single 
multiconductor  electrode  streamer  cable.  Components 
of  this  streamer  include  transmitter  and  receiver  cables, 
connectors,  electrodes,  a  strength  member,  an  outer 
jacket,  and  (in  some  cases)  floats  or  other  buoyancy 
control  means. 

The  gauge  of  the  transmitter  cables  depends  on  the 
maximum  transmitted  current.  Eighteen-gauge  cable  is 
adequate  for  current  levels  of  about  10  A.  Because  the 
receiver  cables  carry  vinually  no  current,  they  can  be  of 
very  light  gauge.  As  mentioned,  current  electrodes  are 
wrapped  around  the  outer  cable  jacket,  and  voltage 
electrodes  are  spliced  to  the  outside  of  the  jacket. 

Integrity  of  cable  and  connector  insulation  is  crucial 
to  data  quality.  Leakage  of  current  to  the  seawater  at 


locations  other  than  the  current  electrodes  can  result  in 
erroneous  voltage  readings;  contact  between  seawater 
and  the  metallic  conductors  leading  to  the  voltage 
electrodes  will  generate  spurious  voltages,  noise,  and 
drift.  Thus  great  care  must  be  taken  to  ensure  that  all 
connections  are  watertight  and  that  no  leaks  exist  in  the 
cable  insulation. 

To  protect  the  individual  cables  from  damage,  they 
are  enclosed  in  an  outer  jacket.  Jackets  of  polyurethane 
material  are  used  for  seismic  cables  and  have  proven 
sufficiently  rugged  for  bottom-towed  operations  in 
which  the  cable  is  dragged  across  the  sea  floor.  A 
central  strength  member  of  Kevlar  or  similar  rope  is 
used  to  relieve  the  mechanical  towing  stress  on  the 
electrical  cables. 

Surface-towed  cables  are  fabricated  with  a  buoyant 
outer  jacket,  and  bottom-towed  cables  are  made 
negatively  buoyant.  A  streamer  cable  to  be  towed  within 
the  water  column  would  be  neutrally  buoyant,  with 
towing  depth  controlled  by  a  V-fin  or  similar  device. 

Instrumentation 

Instrumentation  consists  of  electronic  devices  to 
measure,  record,  and  store  the  transmitted  current  and 
received  voltage  values.  The  on-board  instrumentation 
also  may  perform  initial  data  reduction  and  interpreta¬ 
tion. 

The  simplest  system  for  acquiring  marine  DC  data 
consists  of  a  multichannel  strip  chart  recorder  with 
appropriate  preamplifiers  for  the  individual  current  and 
voltage  channels.  Such  a  system  is  relatively  inexpensive 
and  reliable,  and  allows  the  interpreter  to  directly 
observe  data  quality.  However,  because  a  full  data  set  is 
acquired  with  each  current  pulse  every  few  seconds  of 
operation,  a  great  deal  of  time  is  required  to  later 
manually  read  and  record  the  current  and  voltage  levels 
for  each  transir.ia^d  current  pulse. 

A  computer-controlled  multichannel  digital  data 
acquisition  system  can  be  used  to  greatly  reduce  data 
reduction  and  interpretation  time.  In  such  a  system, 
each  current  and  voltage  pulse  is  automatically  digitized 
and  stored  on  tape  or  floppy  disks  for  later  reduction 
and  interpretation.  This  type  of  system  wruld  be 
considerably  more  complex  and  expensive  than  the  strip 
chart  recorder  described  above,  but  would  greatly 
reduce  dau  interpreuticn  time.  Further  time  savings 
could  be  effected  by  feeding  necessary  supplementary 
information  such  as  vessel  position  and  water  depth  and 
resistivity  directly  into  the  data  acquisition  system. 

Data  reduction  and 
interpretation 

A  marine  DC  resistivity  data  set  consists  of  the 
peak-to-peak  amplitudes  of  sequential  positive  and 
negative  current  pulses  and  the  peak-to-peak  ampli¬ 
tudes  of  the  associated  voltage  pulses  across  electrode 
pairs  A/iN,.  M2N2,  .  .  .  Mn^n-  For  each  pulse,  the 
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peak-to-peak  amplitude  is  divided  by  two  to  give  the 
true  zero-to-peak  amplitude  without  the  effects  of  any 
DC  offsets.  For  each  MN  electrode  pair,  the  apparent 
resistivity  for  the  Schlumberger  array  is  calculated  from: 

»  ir  (V/I)((MN/2)3  -  (AB/2)2]I2(AB/2)  (4) 

where  apparent  resistivity  (Am), 

V  m  voltage  at  potential  e'ectrodes  (mV), 

I  s  current  at  transmitting  electrodes  (mA), 

MN/2  B  one-half  of  potential  electrode  spacing 
(m),  and 

AB12  B  one-half  of  current  electrode  spacing  (m) 

For  manual  interpretation,  the  apparent  resistivity 
data  are  presented  on  double  logarithmic  plots  of 
apparent  resistivity  (p^)  versus  electrode  spacing 
(MN/2)  (Fig.  2b)  (Griffiths  and  King.  1965).  By 
convention,  the  elect'ode  spacing  parameter  is  taken  to 
be  one-half  the  potential  electrode  spacing  (MN/2). 
The  logarithmic  spacing  of  electrodes  yields  uniformly 
spaced  data  points  on  the  electrode  spacing  axis  (Fig. 
2b). 

Manual  interpretation  of  apparent  resistivity  data  to 
obtain  layer  thicknesses  and  resistivities  involves 
comparing  the  measured  resistivity  curve  with  a  set  of 
theoretical  curves,  superimposing  or  overlaying  the  field 
curve  on  a  similarly  shaped  group  of  theoretical  curves 
until  a  fit  is  obtained.  For  the  terrestrial  environment, 
complete  solutions  of  layer  thicknesses  and  re.<iistivities 
for  structures  of  up  to  three  layers  are  readily  achievable 
from  published  curves  for  the  Schlumberger  array 
(Griffiths  and  King,  1965),  and  computer  programs  are 
available  to  generate  resistivity  curve  sets  for  structures 
with  an  arbitrary  number  of  layers. 

Computerized  data  inversion  programs  are  now 
available  to  automatically  determine  the  best-fit  layer 
parameters  tor  a  given  data  set.  With  sophisticated 
programming  and  relatively  consistent  layered  structure, 
the  complete  data  analysis  and  structure  evaluation  can 
be  done  by  computer.  To  date,  when  applying  these 
procedures  for  geotechnical  pu'poses,  the  number  of 
layers  and  their  approximate  thicknesses  and  resistivities 
have  been  estimated  first  by  an  experienced  interpreter 
from  an  examination  of  the  apparent  resistivity  curve. 
These  first  estimates  are  then  used  as  input  to  a 
computer  inversion  program  tu  iteratively  refine  the 
estimates  of  tlie  layer  parameters. 

Data  taken  at  the  sea  surface  can  b~  'nterpreted  using 
one  of  a  number  of  published  or  com  nercially  available 
inversion  programs  written  for  onshore  DC  resistivity 
data  interpretation.  For  an  electrode  array  deployed 
within  the  water  column  or  on  the  sea  floor,  Harding 
Lawson  Associates  has  written  a  forward  modeling 
program  (SUBSOND)  and  an  inversion  program 


(SUBVERT)  based  on  DC  resistivity  surface  modeling 
and  inversion  programs  (SOND  and  INVERT)  written  at 
the  University  of  California,  Berkeley,  modified  using 
expressions  for  buried  electrodes  given  by  Alfano 
(1962).  SUBVERT  will  handle  a  sea-floor  structure  of 
up  to  three  layers  (two  layers  overlying  an  infinitely  deep 
basement),  llie  program  can  be  run  on  a  mainframe 
computer  or  on  an  IBM-PC  or  equivalent.  Inversion  of 
one  data  set  requires  a  few  seconds  or  less  on  a 
mainframe  computer  or  a  few  minutes  on  an 
IBM-PC- AT.  For  production  use,  this  program  could 
be  upgraded  to  handle  more  layers  and  reduce  running 
time. 

Because  acquisition  of  a  data  set  requires  only  a  few 
seconds  but  the  inversion  process  for  the  data  set 
requires  a  few  minutes,  interpreted  data  cannot  be 
displayed  as  rapidly  as  they  are  acquired.  However,  data 
quality  can  be  monitored  by  observation  of  the  apparent 
resistivity  values,  which  can  be  calculated  and  displayed 
immediately  using  Equation  4.  After  data  acquisition  has 
been  completed,  the  interpreter  can  decide  how  many 
and  which  data  sets  are  to  be  interpreted.  Tliis  decision 
is  based  on  the  desired  density  of  data  coverage.  For 
example,  for  a  boat  speed  of  5  knots  and  a  transmitter 
frequency  of  1  Hz,  a  full  data  set  is  acquired  for  each 
meter  of  survey  line.  For  a  desired  lateral  resolution  of 
10  m,  only  every  tenth  data  set  would  have  to  be 
interpreted. 

The  accuracy  of  the  inversion  process  can  be  greatly 
improved  if  one  or  more  of  the  layer  parameters  is 
known  in  advance.  For  marine  measurements,  it  is  easy 
to  independently  determine  the  water  depth  and 
resistivity  at  each  station.  Fixing  these  parameters  in  the 
input  to  the  inversion  program  allows  improved 
resolution  of  the  desired  subseafloor  layer  parameters. 
Thus  continuous  measurement  of  seawater  depth  and 
resistivity  should  be  carried  out  as  part  of  any  marine 
resistivity  survey. 

Signal  and  noise  levels  for 
marine  DC  array 

As  part  of  the  feasibility  assessment  of  using  DC 
resistivity  techniques  for  classifying  marine  sediments  in 
situ,  measurements  were  conducted  in  San  Francisco 
Bay  in  water  depths  of  3-4  m  with  a  towed 
Schlumberger  array  to  directly  measure  noise  levels. 

The  resolution  of  sea-floor  layering  attainable  by  DC 
resistivity  measurements  depends  on  distinguishing 
among  apparent  resistivity  curves  for  various  layered 
structures.  This  ability  is  a  function  of  data  signal-to- 
noise  (S/N)  levels.  That  is,  if  error  bars  on  individual  pg 
field  readings  are  large  relative  to  tire  separation 
between  apparent  resistivity  curves  that  represent 
different  layered  strucUirss,  then  the  interpreter  or 
computer  algorithm  will  not  be  abl  to  distinguish  which 
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of  the  layered  structures  represents  the  true  geology. 
Therefore,  eiqMcted  signal  and  noise  levels  for  actual 
field  measurements  must  be  determined  before  geologic 
resolution  capabilities  of  marine  DC  resistivity  systems 
can  be  determined. 

Noise  levels  depend  on  electrode  properties  and 
operating  conditions.  As  discussed  above,  intrinsic  noise 
levels  for  sutionary  Ag-AgCl  electrodes  separated  by  a 
few  meters  and  immersed  in  calm  seawater  are  a  few 
microvolts.  Thus  under  these  conditions  signal  levels 
greater  than  about  0.1  mV  will  give  S/N  ratios  of  100  or 
more,  penratting  resolutions  of  values  better  than 
1%. 

Noise  levels  generated  by  continuous-towed  opera¬ 
tions  will  depend  on  sea  state,  wave  activity,  bottom 
conditions,  towing  speed,  currents,  and  seawater 
temperature  and  salinity  variations  (Corwin  and  Conti, 
1972;  Corwin,  1976).  Measurements  conducted  in  San 
Francisco  Bay  were  used  to  establish  preliminary  S/N 
levels  for  use  in  this  study.  A  streamer  cable  with 
current  electrode  separation  {ABU)  of  0.25  m  and 
voltage  electrode  separations  (MNI2)  of  1,  20,  and  100 
m  was  towed  at  varying  speeds  on  the  surface,  within  the 
water  column,  and  along  the  floor  of  the  bay  (Fig.  3). 
Note  that  for  construction  purposes  the  100-m  MN 
electrode  pair  is  offset  from  the  center  of  the  array. 
Thus  the  noise  levels  on  this  electrode  pair  will  be 
representative,  but  signal  levels  will  be  12.7  times 
greater  than  those  for  a  concentric  array  with  the  same 
AfA’  spacing. 

Water  depth  in  the  tes,  area  was  about  3-4  m,  and 
measured  water  paramet  ;rs  are  given  in  Table  2.  Wave 
height  varied  from  about  30  to  60  cm  throughout  the 
test  period.  Currents  of  up  to  several  knots  velocity  were 
encountered  along  the  test  profile. 

The  array  was  operated  in  the  following  modes:  (1) 
towed  at  the  sea  surface  at  speeds  of  about  4-S  knots 
and  7-8  knots;  (2)  towed  at  about  2-m  water  depth  at 
the  same  speeds;  (3)  deployed  stationary  on  the  sea 
floor;  and  (4)  towed  along  the  sea  floor  at  about  4 
knots.  In  each  mode,  noise  levels  across  each  of  the 
three  MN  pairs  were  measured  with  no  current 
transmitted  to  allow  maximum  receiver  sensitivity. 
Current  in  the  form  of  5-  or  10-A  square  wave  pulses 
then  was  transmitted  across  the  AB  electrodes  to  allow 
determination  of  S/N  levels. 

The  intrinsic  noise  level  of  a  pair  of  Ag-AgCl 
electrodes  with  MN/2  of  1  m  and  resting  stationary  on 
the  sea  floor  was  found  to  be  of  the  order  of  0.02  mV 
(Fig.  4d).  As  MN/2  increases  to  20  m,  the 
corresponding  noise  level  increased  to  about  0.6  mV 
peak-to-peak  (Fig.  Sb);  and  stationary  noise  levels  at 
lOO-m  MN/2  were  about  2  mV  peak-to-peak  (Fig. 
6d).  Thus  the  stationary  noise  level  is  roughly 
proportional  to  electrode  separation  MN/2. 

When  towed  along  the  surface  at  about  5  knots,  noise 
levels  were  about  0.04  1:jV  {  'r  MN/2  of  1  m,  about  0.8 


Table  2.  Water  parameters  for  DC  array  noise  tests, 
San  Francisco  Bay. 


Depth  (m) 

I'arripsrature 

CO 

Salinity  (ppt) 

Resistivity 

(Dm) 

0.5 

20.0 

19.0 

0.36 

2.0 

IS. 5 

20.5 

0.33 

3.5 

19.0 

14.5 

0.48 

mV  for  MN/2  of  20  m,  and  about  2  mV  for  MN/2  of 
100  m  (Fig.  4a,  Sa,  and  6b).  I  hus  noise  increases  only 
slightly  with  the  array  in  motion.  Increasing  towing 
speed  to  about  8  knots  appeared  to  increase  the 
frequency,  but  not  the  amplitude,  of  the  noise  level. 
Arrays  towed  along  the  sea  floor  showed  noise  levels 
very  similar  to  those  for  surface-towed  arrays. 

Current  vtraveforms  for  the  surface-towed  array  (Fig. 
7)  show  that  current  amplitude  is  resolvable  within 
about  1-3%  from  the  chart  recorder  records.  For  a 
potential  electrode  spacing  (MN/2)  of  1  m  and  applied 
current  (/)  of  S  A,  the  received  signal  level  was  about 
175  mV  (Fig.  4c),  which  when  combined  with  noise 
levels  of  0.04  mV  gives  an  S/N  ratio  of  about  5000  and 
a  corresponding  resolution  of  better  than  1%  for  p^.  For 
MN/2  K  20  m  (Fig.  5d),  a  transmitted  current  of  10  A 
gives  a  signal  level  of  about  2.5  mV  and  a  corresponding 
S/N  ratio  of  about  4  in  the  presence  of  a  large  (0.6  mV) 
noise  pulse  and  about  8  for  more  typical  noise  levels  of 
about  0.3  mV  peak-to-peak.  Increasing  current  levels 
to  at  least  20  A  (as  would  be  done  for  an  operational 

array)  would  given  S/N  ratios  of  about  8  to  16,  with 
corresponding  resolution  of  about  6  to  12%. 

These  results  for  20-m  MN/2  are  useful,  as  this 
spacing  would  be  the  maximum  electrode  spacing  for  a 
sea-floor  towed  array  with  a  subseafloor  depth  of 
investigation  of  about  6  m.  Resolution  of  the  noisiest 
value  for  this  array,  then,  would  be  no  worse  than  about 
12%,  and  resolution  would  improve  rapidly  as  MN/2 


VolUge  (MN)  electrodes:  Sensorex  Model  S22SR  Ap-AgCI 
Current  etoctrodes:  18.flauae  copper  bus  wire  wrapped  around 

cable  (8-mm  diameter),  B^m  length 
Chart  recorders:  Linear  nndel  142,  battery  operated 

Transmitter  power  supply:  3.5-kW  motor  generator,  rectifier<ontrol  unit 
(lO-A  maximum),  manual  switching 

Effective  current 

eloctrode  resistance.  Approx.  50  ohms 


Figure  3.  Configuration  of  test  cable  for  NORDA 
marine  direct  current  resistivity  study. 
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M  SURFACE  TOWED. 
5  kt,  I  -  0 


(b)  SURFACE  VOWED, 
e  kl.  I  -  0 


(C)  SURFACE  TOWED, 
5  kt,  I  *  SA 


«d)  STATIONARY  ON  SEA  FLOOR, 
1-0 


(t)  STATIONARY  ON  SEA  aOOR,  (I)  TOWED  ALONG 

I-5A  SEA  FLOOR.  4  kt,  I  =  0 


(g)  TOWED  ALONG 

SEA  FLOOR.  4  kt,  I  =5A 


Figure  4.  Results  of  noise  tests  in  San  Francisco  Bay,  MN/2  =  1  m.  NORDA  marine  direct  current  resistivity  study. 


Figure  5.  Results  of  noise  tests  in  San  Francisco  Bay,  MN/2  *  20  m,  NORDA  marine  direct  current  resistivity  study. 
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SURFACE  TOWED.  5  Kt 


(c)  SURFACE  TOWED. 

8  kt.  I  s  0 

.J. 

_10  mV 

6  sec 
M 

u _ I _ I _ I _ I _ I _ I _ I _ I 


(e)  STATIONARY  ON  (I)  TOWED  ALONG  SEA  FLOOR, 


Figure  6.  Results  of  noise  tests  in  San  Francisco  Bay,  MN/2  -  103  m,  NORDA  marine  direct  current  resistivity 
study. 
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(b)  SURFACE  TOWED,  C  kt 
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Figure  7.  Transmitted  current  in  noise  tests  in  San  Francisco  Bay,  AB/2  =  0.25  m,  NORDA  marine  direct  cur¬ 
rent  resistivity  study. 
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decreases,  signal  levek  increase,  and  noise  levels 
decrease.  Thus  a  practical  array  with  maximum  MNI2  of 
about  20  m  dragged  along  the  sea  floor  should  give 

#  acceptable  resolution  of  subseafloor  layer  parameters  to 
depths  of  6  m. 

Although  the  received  signal  is  detecuble  for  the 
lOO-m  MNI2  pair,  its  amplitude  of  about  1  mV  is  poorly 
resolved  (Fig.  6a).  As  mentioned  previously,  a  practical 
array  with  MNI2  vahws  greater  than  about  20  m  would 

#  use  a  second  pair  of  current  transmitting  electrodes  with 
considerably  greater  separation  ABi2  to  provide  a  larger 
signal  and  acceptable  resolution  for  large  receiving 
electrode  separations. 

In  this  brief  study,  no  attempt  was  made  to  analyse 
the  sources  of  electrode  noise  (which  probably  include 

0  wavv  and  current  activity,  changes  in  seawater 

properties,  and  motion  of  the  connecting  cables  through 
the  earth’s  magnetic  field).  Better  undersunding  of 
these  sources  could  produce  methods  to  reduce  this 
noise  and  thus  improve  S/N  ratios  without  the  necessity 
for  larger  transmitter  currents. 

Analytical  studies  of  DC  array 
performance 

Description 

#  A  series  of  computer  modeling  studies  were 
conduaed  using  the  program  SUBSOND  to  obtain  a 
preliminary  idea  of  the  resolution  of  sea-floor  sediment 
properties  provided  by  tire  marine  DC  resistivity 
method.  The  sea  floor  was  modeled  a.«  a  two-layer 
system:  a  surficial  sediment  layer  over  a  more  resistive 

#  layer  of  infinite  depth.  In  this  study,  the  resolution  of 
the  resistivity  of  the  upper  layer  was  the  parameter  of 
major  interest,  the  thickness  of  the  upper  layer  was  of 
secondary  interest,  and  resistivity  of  the  underlying 
strata  was  considered  the  least  important. 

The  nomenclature  and  range  of  parameters  used  in 

^  these  preliminary  parameter  studies  are  illustrated  in 
Figure  8  and  are  listed  as  follows: 

d  «  water  depth  (10,  100,  and  1000  m  used). 

h  B  height  of  the  electrode  array  off  the  sea  floor  (0. 
1,  10,  30,  and  100  m  above  the  sea  floor  all  used), 

« thickness  of  the  upper  layer  of  the  sea  floor  (1,6, 
60  m  used), 

«  thickness  of  the  second  layer  (infinite). 

Pm.  b  resistivity  of  the  overlying  water  (assumed  0.25 
flm), 

pt  «  resistivity  of  the  upper  layer  of  sea  floor  (0.4, 
0.8,  1.6  Dm  U3ed),  and 

^  Pa  *  resistivity  of  the  second  sea-floor  layer  (2.0  dm 

used). 

The  maximum  potential  electrode  half-separation 
(MSI2)  used  for  all  cases  was  316  m.  Table  3 


DC  ELECTRODE  ARRAY 

_ _ _ _ 

k 

9m  ■  0.2S  ohm-m 
d  -  10, 100, 1000  m 

i _ 

‘h  SEA  FLOOR 

t!] _ - ^  ' 

>  1,  6.  60m  Q■^  m  0.4,  0.8, 1.6  ohm-m 

*  gj  ■  2.0  ohm-m 

Figure  8.  Nomenclature  used  in  de.’.cribing  case 
geometries  for  analytical  studies,  marine  DC  resistivity 
study. 


summarizes  parameters  selected  for  each  of  the  cases 
examined. 

The  ability  to  differentiate  between  a  given  pair  of 
resistivity  curves  depends  on  the  S/N  ratio  of  the 
individual  data  points  defining  the  curves.  The  S/N  ratio 
will  vary  depending  on  operating  conditions,  transmitted 
current  amplitude,  and  voltage  receiver  electrode 
separation  (or  half-spacing,  MNI2).  Thus  no  single  S/N 
ratio  can  be  assumed  for  the  entire  range  of  a  given 
resistivity  curve,  nor  can  the  S/N  ratio  be  assumed  to  be 
the  same  for  different  environments.  Because  this  study 
is  ortly  a  preliminary  effort,  the  S/N  ratio  was  assumed 
to  be  a  constant- IS -for  all  receiver  electrode  spacings, 
or.  said  differently,  a  resolution  level  of  *7%  ha:  been 
assumed  for  all  receiver  electrode  spacings.  Quite 
conveniently,  the  line  width  on  the  computer  plots 
generated  for  this  study  Just  happer*d  to  also  represent 
a  resolution  of  about  7%;  thus,  'k  len  two  computer- 
generated  curves  are  visually  distinguishable,  then 
field-measured  data  should  also  be  resolvable  to  the 
accuracies  required.  (For  example.  Figure  A-1  of 
Appendix  A  presents  the  computer-generated  perform¬ 
ance  prediction  for  Case  1  (see  Table  3),  showing  that 
Layer  1  resistivities  are  resolvable  for  the  1-m-thick 
layer  only  for  receiver  electrode  half-spacing  (MNI2) 
less  than  3  m. 

S/N  ratios  will  be  considerably  greater  than  IS  for 
most  of  the  inner  potential  elecuodes  (MN)  under  most 
opsrating  conditions;  therefore,  computer-generated 
curves  that  were  not  visually  distinguishable  at  small 
receiver  spacings  {MNI2)  often  may  be  resolvable  from 
field  data.  Therefore,  using  15  for  the  value  of  the  S/N 
ratio  over  the  e’ltire  calculated  curve  is  a  conservative 
assumption. 

Results 

The  computer-generated  plots  of  apparent  resistivity 
(Pa)  versus  teceiver  electrode  half-spacing  {MN/2)  for 
each  of  these  cases  listed  in  Table  3  are  reproduced  in 
Appendix  A  (from  Corwin,  1986).  The  essential 
contents  of  the  Appendix  A  figures  are  summarized  in 
Table  4,  where  the  electrode  spacing  (MN/2)  required 
to  resolve  the  resistivity  of  the  sediment  layer  1  (p^)  is 
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Tabit  3.  Layer  parameters  for  model  studies  (resolution  of  sea-floor  resistivity). 


(Om) 


Cate 

a 

(m) 

h 

(ml 

‘I 

(ml 

Curve 

1 

Curve 

2 

Curve 

3 

P2 

(flmi 

1 

10 

0 

1 

0.4 

0.8 

1.6 

2.0 

2 

10 

0 

6 

0.1 

0.8 

1.6 

2.0 

3 

10 

0 

60 

0.4 

0.8 

1.6 

2.0 

4 

10 

1 

1 

0.4 

0.8 

1.6 

2.0 

5 

10 

1 

6 

0.4 

0.8 

1.6 

2.0 

6 

10 

1 

60 

0.4 

0.8 

1.6 

2.0 

7 

10 

10 

1 

0.4 

0.8 

1.6 

2.0 

8 

10 

10 

6 

0.4 

0.8 

1.6 

2.0 

8 

10 

10 

60 

0.4 

0.8 

1.6 

2.0 

10 

100 

0 

1 

0.4 

0.8 

1.6 

2.0 

11 

100 

0 

6 

0.4 

0.8 

1.6 

2.0 

12 

100 

0 

60 

0.4 

0.8 

1.6 

2.0 

13 

100 

1 

1 

0.4 

0.8 

1.6 

2.0 

14 

100 

1 

6 

0.4 

0.8 

1.6 

2.0 

15 

ICO 

1 

60 

0.4 

0.8 

1.6 

2.0 

16 

100 

30 

1 

0.4 

0.6 

1.6 

2.0 

17 

100 

30 

6 

0.4 

0.8 

1.6 

2.0 

18 

100 

3C 

60 

0.4 

0.8 

1.6 

2.0 

19 

100 

100 

1 

0.4 

0.8 

1.6 

2.0 

20 

100 

100 

6 

0.4 

0.8 

1.6 

2.0 

21 

100 

100 

60 

0.4 

0.6 

1.6 

2.0 

22 

1000 

0 

1 

0.4 

0.6 

1.6 

2.0 

23 

1000 

0 

6 

0.4 

0.8 

1.6 

2.0 

24 

1000 

0 

60 

0.4 

0.8 

1.6 

2.0 

25 

1000 

1 

1 

0.4 

0.8 

1.6 

2.0 

26 

1000 

1 

6 

0.4 

0.8 

1.6 

2.0 

27 

1000 

1 

60 

0.4 

0.8 

1.6 

2.0 

28 

1000 

30 

1 

0.4 

0.8 

1.6 

2.0 

29 

1000 

30 

6 

0.4 

0.8 

1.6 

2.0 

30 

1000 

30 

60 

0.4 

0.6 

1.6 

2.0 

31 

1000 

100 

1 

0.4 

0.6 

1.6 

2.0 

32 

1000 

100 

6 

0.4 

0.8 

1.6 

2.0 

33 

1000 

100 

60 

0.4 

0.8 

1.6 

2.0 

listed  for  each  case.  The  data  of  Table  4  is  rearranged  in 
Table  S,  flrst  in  terms  of  the  thickness  of  the  surficial 
sediment  layer,  second  in  terms  of  the  height  of  the 
array  above  the  sea  floor,  and  third  in  terms  of  the 
water  depth  (to  the  sea  floor). 

Table  S  shows  that  the  resistivity  of  a  1-m-thick 
surficial  sediment  layer  can  be  resolved  as  being  0.4, 
0.8,  or  1.6  Om  only  when  the  array  is  laying  on  the  sea 
floor  and  oitly  for  receiver  electrode  half-spacings 
{MNI2)  of  3  m  or  less.  Even  when  the  array  is  towed  as 
close  as  1  m  above  the  sea  floor,  the  resolution  of  the 
apparent  resistivity  (pg)  for  the  surficial  1-m  layer  is  not 
adequate  to  permit  individual  sediment  classification  for 
that  1-m  layer. 

When  the  thickness  of  the  surficial  sediment  layer  is  6 
m,  then  sufficient  resolution  can  be  achieved  to  classify 


tlie  sediments  of  that  layer  by  towing  the  array  up  to  10 
m  above  the  sea  floor.  When  towed  at  10  m  above  the 
sea  floor,  the  half-spacing  of  one  or  more  pairs  of 
receiving  electrodes  must  be  at  least  12  m,  requiring  an 
array  length  of  at  least  24  m,  in  order  to  adequately 
resolve  the  apparent  resistivity  (p^)  and  to  permit 
calculation  of  the  layer  1  resistivity  (pO- 

For  a  surficial  sediment  layer  60  m  thick, 
classification  is  possible  with  an  array  towed  anywhere 
from  on  the  sea  floor  to  30  m  above  it.  However,  as  the 
position  of  the  array  is  raised  above  the  sea  floor,  the 
spacing  of  receiver  electrodes  required  to  achieve  the 
required  resolution  for  sediment  classification  also 
increases.  For  an  array  30  m  above  the  sea  floor  the 
required  electrode  half-spacing  varies  with  the  water 
depth  (</):  for  d  of  100  m,  MN/2  must  be  at  least  35  m, 
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Table  4.  Electrode  spacInQs  (MN/2)  required 
in  order  to  resolve  pv 


a 

Cast  (m) 

h 

(m) 

(m) 

MN/2  Required 
(m) 

1 

10 

0 

1 

^3 

2 

10 

0 

6 

^140 

3 

10 

0 

60 

1  through  316 

4 

10 

1 

1 

None 

5 

10 

1 

6 

2  through  130 

e 

10 

1 

60 

1  through  316 

7 

10 

10 

1 

None 

e 

10 

10 

6 

12  through  130 

9 

10 

10 

60 

6  through  316 

10 

100 

0 

1 

1  through3 

11 

100 

0 

6 

1  throughlS 

12 

100 

0 

60 

1  through  316 

13 

100 

1 

1 

None 

14 

100 

1 

6 

1  through  22 

1$ 

too 

1 

60 

1  through  316 

16 

100 

30 

1 

None 

17 

100 

30 

6 

None 

16 

100 

30 

60 

35  through  316 

18 

100 

100 

1 

None 

20 

100 

100 

6 

None 

21 

100 

100 

60 

110  through  316 

22 

1000 

0 

1 

1  through  30 

23 

1000 

0 

6 

1  through  15 

24 

1000 

0 

60 

1  through  ISO 

26 

1000 

1 

1 

None 

26 

1000 

1 

6 

1  through  20 

27 

1000 

1 

60 

1  through  ISO 

26 

1000 

30 

1 

None 

28 

1000 

30 

6 

None 

30 

1000 

30 

60 

60  through  200 

31 

1000 

100 

1 

None 

32 

1000 

100 

6 

None 

33 

1000 

100 

60 

None 

and  for  d  of  1000  m.  MN/2  must  be  60  m.  It  is 
significant  to  note  that  an  array  towed  100  m  above  the 
sea  floor  achieves  sufficient  resolution  to  distinguish 
resistivities  (  for  a  60-m-thick  surficial  layer  when 
the  water  depth  is  100  m,  but  does  not  have  the 
required  resolution  in  1000-m  water  depth. 

Discussion 

The  formation  factors  u'^d  for  these  cases  (1.6,  3.2. 
and  6.4)  (resistivities  of  0.4,  0.8,  and  1.6  dm)  typically 
would  represent,  respectively,  a  poorly  consolidated, 
fine'grained  sediment  (e.g.,  clay  or  fine  silt);  a 
moderately  consolidated,  medium-grained  material 
(e.g.,  fine  sand  or  silty  sand);  and  coarse-grained, 
well-compacted  material  (e.g.,  gravel  or  compacted 
sand).  Thus  the  rcsolvability  of  the  three  curves  on  each 


figure  of  ^}pendix  A  is  indicative  of  our  ability  to 
distinguish  iMtween  the  most  common  sea-floor 
materials. 

Con^rison  between  cases  with  values  of  1,  6,  and 
60  m  indicates  that  it  is  increasingly  easy  to  resolve  the 
resistivity  of  the  upper  layer  as  its  thickness  increases. 
For  all  of  the  cases,  in  any  water  depth,  an  array  on  or 
just  above  the  sea  floor  can  resolve  the  resistivity  of  a 
1-m-thick  tqpper  layer.  An  array  at  the  sea  surface  in 
10  m  of  water  cannot  resolve  the  resistivity  of  a 
1-m-thick  layer  (Fig.  A-7,  Table  5),  but  can  resolve 
the  resistivity  of  a  6-  or  60-m-thick  upper  layer  (Figs. 
A-8  and  A-9,  Table  5). 

For  100  m  or  1000  m  water  depth,  an  array  30  m 
above  the  tea  floor  or  at  the  surface  carmot  resolve  the 
the  resistivity  of  a  1-  or  6-m-thick  upper  layer  but  can 
resolve  tliat  of  a  60-m-thick  layer.  Thus  for  any  water 
depth  the  array  must  be  towed  at  or  close  to  the  sea 
floor  to  resolve  the  resistivity  of  layers  that  are  thin 
relative  to  the  water  depth.  As  the  thickness  of  a  given 
layer  becomes  larger  relative  to  water  depth,  its 
resistivity  becomes  easier  to  determine,  and  the  array 
can  be  towed  at  a  greater  height  off  the  sea  floor. 

Electromagnetic  methods 

The  technique 

The  EM  method  proposed  for  measuring  the 
electrical  resistivity  of  surficial  sea-floor  sediments  is 
(Hie  of  a  family  of  geoelectrical  prospecting  techniques 
that  utilize  a  controlled  time-varying  current  to  generate 
a  source  field.  These  primary  fields  induce  currents 
irithin  the  surrounding  media.  The  amplitude  and  phase 
of  these  secondary  fields  reflect  the  variation  in  physical 
parameters  that  occurs  in  the  vicinity  of  the  current 
source.  This  type  of  eiqiloration  method  will  be  referred 
to  as  the  EM  teennique.  The  EM  methods  described 
herein  use  a  current  source,  which  varies  over  a  band  of 
frequencies  from  10  to  10,000  Hz. 

In  implementing  the  method,  a  magnetic  field  is 
generated  by  driving  a  time-varying  electric  current 
through  a  loop  of  wire  or  through  a  straight  wire 
grounded  at  both  ends.  When  any  conductive  materials, 
such  as  seawater  and  sea-floor  sedimtnts,  are  present, 
induced  eddy  currents  will  flow  in  closed  paths  or  loops 
normal  to  the  direction  of  magnetic  field.  These  eddy 
currents,  in  turn,  generate  secondary  magnetic  fields. 
Thus,  at  any  point  in  qiace,  the  total  magnetic  field  may 
be  thought  of  as  consisting  of  two  parts:  (1)  a  primary  or 
normal  field  due  to  the  source  current  and  (2)  a 
secondary  or  disturbing  field  due  to  eddy  currents 
induced  in  conductors  (Keller  and  Frischknecht,  1966). 

The  functioning  of  an  EM/induction  system  may  be 
compared  with  that  of  •  uansformer  having  three 
windings,  with  one  of  the  windings  being  shorted  to 
simulate  the  effect  of  a  conductive  body  (see  Fig.  9).  A 
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Table  6.  Paramatart  of  the  seven  cases  examined  In  the  EM  analytical  studies. 


Ca«e 

d 

(m) 

h 

(m) 

t, 

|m) 

p« 

(Oin) 

-P* 

(flm) 

1A 

oo 

10 

oo 

2.50 

IB 

oo 

10 

oo 

'..2S 

— 

1C 

00 

10 

oo 

0.60 

2A 

oo 

10 

e 

1.25 

2.5 

2B 

oo 

10 

e 

0.60 

2.5 

SA 

oo 

11 

oo 

2.50 

SB 

oo 

9 

oo 

2.50 

4A 

15 

10 

oo 

2.50 

— 

SA 

20 

10 

oo 

2.50 

SB 

20 

10 

oo 

1.25 

— 

SC 

20 

10 

oo 

0.60 

... 

SA 

oo 

OO-10 

- 

— 

— 

7A 

oo 

20 

oo 

2.50 

7B 

oo 

30 

oo 

2.50 

... 

soiid  conductive  body  (such  as  the  seawater  column  or 
the  ledimcnts  beneath)  behaves  in  a  simila'  way  to  the 
shorted  turn,  except  that  the  distribution  of  eddy 
currents  depends  not  ortly  on  frequency  but  also  on  the 
shape,  size  (thickness  t  for  horizonul  layers),  and 
electrical  resistivity  (p^)  of  the  body.  At  low  applied 
frequencies  eddy  currents  flow  deep  within  the 
conductor.  At  high  frequencies,  the  skin  effect  restricts 
current  flow  to  a  region  near  the  surface  of  the  body. 

EM  fields  depend  on  the  magnetic  permeability  of 
earth  materials,  as  well  as  on  Ute  electrical  conductivity. 
Ordinarily  the  range  of  effects  that  can  be  attributed  to 
variations  in  magnetic  permeability  is  small  compared  to 
the  range  of  effects  caused  by  variations  in  electrical 
conductivity.  However,  nearby,  highly  magnetic  rocks 
may  distort  this  normal  situation,  altering  the  usual 
relationship  between  primary  and  secondary  fields. 
These  perturbations  are  normally  small  in  a  marine 
environment,  but  a  variable  permeability  may  be 
incorporated  into  the  interpretation  procedure  when 
required. 

The  marine  envirorunent  offers  Kveral  advanuges  to 
the  use  of  EM  techrtiques  that  are  not  enjoyed  on  land. 
The  highly  conductive  seawater  piuvides  a  low 
impedance  contact  for  the  transmitter  antenna  with  no 
more  hardware  involved  than  a  bare  wire  in  the  water. 
Because  of  the  low  impedance,  a  very  large  transmitter 


Figure  9.  Tran^ormer  analogy  of  an  induction  system, 
a  primary  magnetic  field  Hp  is  generated  by  an 
oscillatory  current  flowing  in  the  source-coil.  This  field 
generates  current  in  a  shorted  turn,  representing  a  con¬ 
ductive  earth.  This  current  generates  a  secondary  magnetic 
field.  H^.  Both  fields  are  measured  by  the  voltage 
induced  in  a  receiver  coil.  (From  Keller  and  Frischknecht, 
1966) 
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current  of  30C  A  may  be  obuinad  from  a  moderately 
powerful  aupply  (Corttubic  et  al.,  1986).  Dipole-type 
anteruut  may  be  towed  through  the  eater  without 
interrupting  the  electrical  conuct.  For  a  receiving 
antenna  towed  near  the  ma  floor,  much  of  the 
ionoqjheric  noUe  is  removed  in  the  thickness  of 
overlying  water,  which  makes  the  resolution  of  very 
small  tigitals  possible.  Last,  since  the  seawater  is  more 
conductive  than  the  sea-floor  materials,  at  the 
frequencies  of  operation  (10  Hz  to  1  kHz)  much  of  the 
signal  propagates  through  the  sea-floor  materials  (the 
region  of  interest)  and  not  through  the  seawater  column. 
In  very  shallow  water  where  the  antennas  cannot  be 
positioned  a  large  distance  from  the  air-sea  interface, 
tome  of  the  EM  signal  does  ukc  an  atmospheric  path, 
but  this  path  is  unimportant  when  the  water  depth  is 
more  than  one  skin  depth  (2S  m  at  100  Hz)  (Constable 
et  al..  1986). 

The  major  limitation  on  the  use  of  underwater 
electromagnetic  systems  arises  as  a  resAiIt  of  the  high 
attenuation  of  electromagnetic  fields  as  they  pn^gate 
through  the  seawater  and  the  highly  conductive  surficial 
sediments.  This  attenuation  is  usually  presented  in  terms 
of  a  skin  depth; 


where 
w  «  2trf 

/  >  signal  frequency  (Hz), 
p  «  magnetic  permeability  (H>nry/m),  and 
a  c  electrical  conductivity  (Siemens/m),  o  n  //p 
As  an  exam(.le.  a  unit  amplitude  100-Hz  signal, 
without  considering  geometric  qireading.  woi'ld  be 
attenuated  as  it  propagates  through  the  seawater  (04,  « 
4  S/m)  to  e-i  of  iu  origin  il  amplitude  in  propagating 
approximately  25  m.  (See  Fig.  10).  Thuc.  the  high 
attenuation  of  signal  will  require  that  £M  systems  used 
for  sea-floor  description  sediment  classificaticn  be 
deplo)ed  as  close  as  possible  to  the  sea  flook'. 

Array  geometries 

Straight-wzre  system 

Constable  et  al.  (198i.)  describe  a  straight-w'-e  dipole 
system  with  grounded  enns  intended  for  use  in 
determining  geologic  structure  in  mvrire  areas  where 
seismic  reflection  techniques  perform  poorly,  'ilie 
dipole-dipole  system  developed  consists  of  a  30-m 
transmitting  dipole  towcu  in  contact  with  the  sea  floor, 
injecting  300  A  into  the  seawater  at  frequencies  of  10 
Hz  to  1  kHz  (Figs.  11a  and  b).  Power  is  transmitted 
down  the  towing  cable  at  1  kV  to  limit  losses  due  to 
cable  resisance  (Fic.  lib).  A  transformer  near  the  sea 
floor  convers  the  power  to  low  voltage,  high  current. 


FREQUENCY  (Hr) 


Fifure  10.  Skin  depth  of  electromagnetic  signal  as  a 
function  of  signal  frequency  and  water  conductivity 
(conductivity  is  the  reciprocal  of  resistivity). 


Figure  11.  Towed  EM  system.  ’'A‘*  shows  the  geometry 
of  an  in-line  towing  system.  Black  discs  represent  acoustic 
transponders  used  for  navigation.  "B"  and  "C"  show 
schematics  of  the  transmitter  and  receiver,  respectively. 
(After  Constable  et  al.,  1986) 


Dipole  receivers  with  Ag-AgCl  electrodes  are  (1)  towed 
in  undem  with  the  transmitter  from  the  same  ship.  (2) 
tewed  at  a  disunce  firom  a  second  ship,  or  (3)  deployed 
as  stationary  receivers  on  the  sea  floor.  Data  are 
collected  on  magnvuc  tape  and/or  transmitted  by  radio 
to  the  u  ansmitter  ship  (Fig.  11c).  Positions  of  the  dipole 
arrays  arc  esublished  by  means  of  acoustic  transponders 
fitted  to  the  transmitter,  the  receiver,  and  the  ship  (Fig. 
11). 

The  dipole-dipcle  system  includes  special  receivers 
that  are  capable  of  measuring  very  small  signals  and  that 
take  advantage  of  the  low  input  impedance  offered  by 
the  seawater  (Webb  et  al..  1985).  At  the  highest 
frequencies,  the  receiver  electrode  system  has  intrinsic 
noise  levels  as  low  as  the  thermal  agitation  noise  of  the 
3  n  electrode  grounding  resistance  (4  x  10'*®  V*/Hz 
above  30  Hz).  At  frequencies  below  1  Hz.  electrode 
noise  increases  roughly  proportional  to  the  inverse  of 
frequency.  Receiver  noise  in  the  towed  array  is 


dominated  by  60-Hz  (^iwerline  noise  and  ionospheric 
noise.  Towing  of  the  receiven  at  2  knots  does  result  in 
an  increase  in  noise  at  frequencies  above  GO  Kt,  but  the 
increase  is  not  oveipowering.  At  frequencies  below  60 
Ha,  towing  increases  noise  from  10  to  100  times  due  to 
movement  of  the  receiver  antenna  (Constable  et  at.. 
1986). 

Loop-Loop  Systems 

Loop-Loop  electromagnetic  systems  are  available 
commercially  for  search  and  localisation  for  objects  that 
are  either  lost  on  or  are  shallowly  buried  in  the  sea 
floor.  These  search  systems  are  designed  to  detect  small 
electrical-conducting  (meul)  masses,  and  their  hard* 
ware  design  is  likely  not  optimal  for  the  delineation  of 
horiaontal  boundaries  of  sediment  layers  and  measure¬ 
ment  of  electrical  resistivities  of  those  layers. 

Airborne  electromagnetic  (AEM)  loop-loop  systems 
have  been  used  to  measure  sea  ice  thickness  and  water 
depth  to  the  sea  floor,  and  to  measure  resistivities  of 
surficial  sea-floor  materials.  The  AEM  work  at  NORDA 
has  used  a  commercially  available  Dighem  III 
frequency-domain  system  (described  by  Fraser,  1978, 
1979,  and  1981)  without  significant  modificarion.  The 
Dighem  III  AEM  system  uses  two  horizonul-coplanar 
coil  pairs  sefiarated  by  8  m  and  operating  at  38S  Hz  and 
7200  Hz.  The  transmitter-sensor  (coil  pair)  was 
mounted  in  a  bird  and  towed  by  a  helicopter  at  the  end 
of  a  30-m  cable,  with  the  sensor  maintained  between  40 
and  50  m  above  the  sea  surface.  The  resultant  AEM 

field  dau  from  track  lines  in  the  Cape  Cod  Bay  yield 
reasonable  predictions  of  sea-floor  sediment  conductiv¬ 
ity  up  to  one  skin  depth  of  the  source  frequency:  for  the 


Figure  12.  Coil  cor\figuration5  for  two-loop  soundings. 
(From  Keller  and  Frischknecht,  1966) 


employed  frequencies  of  385  Hz  and  7200  Hz  and  for 
the  given  water  electrical  resistivity  of  0.25  Dm,  the  skin 
depths  are  12.8  m  and  2.9  m.  reqjectively  (Won  and 
Smiu,  1985). 

In  the  design  of  a  loop-loop  system,  the  separation 
between  the  lourcti  loop  and  the  receiver  loop  must  be 
somewhat  more  than  five  times  the  diameter  of  either 
loc^  in  order  to  facilitate  the  arulysis  and  permit  both 
loops  to  be  treated  as  dipoles  (Wait.  1954).  Four  ways 
are  commonly  used  in  which  the  transmitting  and 
receiving  loops  are  oriented  with  respect  to  each  other 
for  deteimining  the  thickness  and  resistivity  of  layered 
material:  horizonul  coplanar.  vertical  coplanar,  vertical 
coaxial,  and  perpen^cular  (Fig.  12)  (K.eller  and 
Frischknecht.  1966). 


Development  of  predictive 
techniques 

Solutions  for  the  vertical  sounding  problem  using  EM 
teclmiques  in  layered  terrestrial  materials  have  been 
developed,  and  computer-generated  solutions  are 
available.  In  addition,  a  number  of  inversion  techniques 
have  been  developed  for  the  interpreution  of  dau 
.'’cquircd  on  the  surface  of  a  layered  earth  (Fullagar  and 
Oldenburg  1984;  Olenn  ct  al.,  1973). 

As  part  of  the  first-year  effort  on  this  project, 
software  has  been  developed  to  provide  accurate 
simulations  of  the  EM  response  for  various  antennas  in 
a  conduciing  media.  The  theoretical  model  used  is 
based  on  the  work  of  Kong  (1972)  and  Tang  (1979). 
This  mathematical  approach  provides  the  closed-form 
solution  of  Maxwell's  equations  for  arbitrary  dipole 
sources  in  an  n-layered  anisotropic  media.  The  solution 
is  decomposed  into  two  modes.  Transverse  Electrical 
(TE)  and  Transverse  Magnetic  (TM).  These  model 
equations  provide  a  set  of  integral  equations,  which  have 
been  solved  using  an  adaptive  numerical  intergration 
package  developed  by  Chave  (1983).  A  portion  of  the 
model  equations  developed  and  validated  as  part  of  this 
effort  are  presented  in  Appendix  B. 

Software  for  a  VAX  11-780  was  developed  from  the 
EM  model  equations  of  Appendix  B  and  was 
successfully  subjected  to  a  number  of  internal  checks.  In 
addition,  the  results  of  the  NORDA  software  have 
checked  well  with  the  published  results  of  Coggen  and 
Morrison  (1970)  for  the  vertical  magnetic  dipole  case 
and  those  of  Weaver  (1967)  for  electric  dipoles  buried 
in  a  simple  two-layer  model.  The  model  analyses  done 
for  this  study  were  calculated  on  the  VAX  1 1-780,  and 
double-precision  accuracy  was  used  for  ai'  numerical 
manipulations.  The  integrations  were  calculated  using  an 
absolute  error  criteria  of  10*^^. 
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Analytical  studies  of  EM  array 
performance 

Description 

Using  the  software  described  above,  electric  and 
magnetic  fields  were  calculated  for  seven  basic  cases. 
These  simula:io>i.  provided  k'  suiting  point  for 
determining  the  rvlative  iibility  of  various  source  types  to 
detect  expected  variations  in  s.^a-floor  sediment 
resistivities.  The  effecu  of  variations  in  source 
antenna-receiver  antenna  roll  and  pitch  also  have  been 
examined.  The  effects  of  the  air-wate.''  interface  on 
resistivities  measured  are  considered,  end  the  interac* 
tion  between  the  secondary  Helds  generated  by  both  the 
water  and  sediment  surfaces  will  be  discussed. 

The  seven  cases  examined  were  designed  to  answer 
several  questions  concerning  antenna  deployment, 
system  characteristics,  and  dau  reduction  requirements. 
A  set  of  sketches  describing  the  various  case  geometries 
is  provided  in  Appendix  C.  In  the  cases  considered,  the 
transmitter  and  receiver  antennas  were  assumed 
atuched  to  a  rigid  boom  with  a  separation  center-to- 
center  oi  10  m.  This  separation  of  10  m  was  dicuted  by 
the  assumption  of  a  practical  diameter  of  1  m  for  tlie 
transmitter.  A  ratio  of  antenna  sepiiration-to-transmit* 
ter  diameter  of  10  or  more  is  assumed  necessary  to 
obuin  valid  results  for  the  analytical  models  developed 
(Appendix  B),  because  these  models  are  derived 
assuming  that  the  transmitter  is  an  idealized  dipole 
source. 

Electric  field  values  for  this  EM  case  study  are  in 
units  of  volts  per  meter  (V/m)  with  a  source  moment  of 
unity.  Magnetic  field  values  are  for  unit  moment  sources 
and  are  in  units  of  Ampere-tums  per  meter 
(A-tums/m).  The  magnetic  permeabilities  (p)  of  all 
media  were  assumed  to  be  x  10-^  Heiuys  per  meter 
(H/m).  and  the  magnetic  field  values  presented  may  be 
converted  to  magnetic  induction  units  by  using 

1  Ampere-tum/meier  ■  1.2566  x  1(P  nT  (y). 

The  electric  and  magnetic  field  values  are  composed 
of  contributions  from  Primary  and  Secondary  fields  in 
the  form  of  Real  (In  Phase)  and  Imaginary  (Quadra¬ 
ture)  Components.  The  model  data  are  presented  in  the 
next  section;  these  electric  and  magnetic  field  values  are 
plotted  on  a  linear  scale  versus  an  induction  number  on 
a  log  scale: 

Induction  number  «  (wpa  )V*r,  (6) 

where 

tt  ■  211/, 

/  >  frequency  (Hz),  range  examined  from  10  to 
10,000  Hz. 

4  «  magnetic  permeability  (H/m),  set  at  4ir  x  10'^ 
H/m, 


a  ■  electrical  conductivity  (S/m),  set  at  4  S/m,  where 
a  «  l/p, 

p  ■  elsictrical  resistivity  (Am),  set  at  0.25  Am,  and 

r  ■  disunce  between  centers  of  transmitter  and 
receiver  (m),  set  at  10  m. 

The  technical  feasibility  of  the  underwater  EM 
sediment  classifleation  concept  will  largely  depend  on 
the  noise  level  present  in  the  measurement  environ¬ 
ment.  For  purposes  of  this  preliminary  study  of  EM 
performance,  noise  level  criteria  have  been  developed 
from  two  sources  of  field  data.  First,  the  electric  field 
noise  data  presented  earlier  (Corwin,  1986),  believed 
valid  for  froqusneies  below  1  Hz,  are  expected  to 
represent  an  upper  limit  for  the  noise  level,  with  electric 
fields  at  higher  frequencies  characterized  by  lower  noise 
levels  due  to  increased  propagation  attenuation  with 
frequency.  Corwin  observed  noise  levels  of  2  to  5  pV 
per  meter  for  his  DC  resistivity  Schlumberger  array; 
thus,  for  the  EM  study  an  electric  field  noise  level  of  5 
pV/m  has  been  assumed.  The  magnetic  field  noise  level 
was  selected  based  on  dau  from  an  experiment 
conducted  by  Oiannini  and  Thayer  (1982).  Giannini 
and  Thayer  towed  a  horizontal  electric  dipole  source  at 
depths  below  the  sea  surface  of  20  to  96  m  at  a  speed  of 
3  m/sec .  and  meesured  magnetic  fields  over  depths 
below  the  surface  from  46  to  109  m.  In  the  frequency 
range  of  1  to  10  Hz,  the  noise  level  was  determined  to 
be  about  0.002  nT.  which  is  equivalent  to  a  magnetic 
field  level  of  1.5  x  10**  A-tum/m.  Based  on  these  data, 
a  magnetic  field  noise  level  of  2  x  10**  A-tum/m  was 
used  in  this  study. 

Results 

System  10  m  above  sea  floor,  deep  water 

Figures  13a  and  b  present  EM  system  performance 
results  for  an  antenna  system  positioned  as  shown  in 
Figure  C-1  (of  Appendix  C)  with  a  transmitter  moment 
of  100  A-tum-m*.  The  antenna  system  is  shown 
operating  10  m  above  the  sea  floor  in  very  deep  water 
over  three  sediment  types  with  resistivities  of  0.62,  1.25 
and  2.5  Am  (i.e..  Formation  Facton  of  2.5,  5,  and  10). 
Calculations  have  been  performed  for  horizontal  and 
vertical  dipole  sources.  The  use  of  a  horizontal  dipole 
source  provides  both  Transverse  Electric  (TE)  and 
Transverse  Magnetic  (TM)  coupling  into  the  sediments 
and  theoretically  provides  more  information  on 
sediment  resistivities  than  the  other  source  types.  The 
vertical  total  magnetic  field  generated  by  a  horizontal 
magnetic  dipole  is  shown  in  Figures  13a  and  b,  with  real 
or  in-phase  response  in  Figure  13a  and  the  imaginary  or 
quadrature  components  given  in  Figure  13b.  The  real 
and  imaginary  responses  have  units  of  A-turn/m  and 
are  plotted  on  an  arithmetic  scale  on  the  ordinate  or 
vertical  axis  of  the  figures.  Tnese  responses  are  plotted 
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Figure  13(a).  Resoiution  of  sediment  types,  array  10  m 
above  seafloor,  overlain  by  a  very  thick  layer  of  water; 
real  component  of  vertical  magnetic  field  (Hz), 
generated  by  horizontal  magnetic  dipole  (HMD). 


Figure  13(b).  Resolution  of  sediment  types,  array  10  m 
above  sea  floor,  overlain  by  a  thick  layer  of  water, 
imaginary  component  of  Hz  generated  by  HMD. 


against  corresponding  values  of  induction  number  (see 
Equation  6)  on  a  logarithmic  scale,  where  the  induction 
number  is  a  function  of  frequency  which  varies  from  10 
to  10,000  Hz. 

Given  expected  equipment  performance,  variations 
on  the  order  2x10*^  A-tum/m  in  magnetic  field  should 
be  developed  by  a  source  of  unit  moment,  i.e.,  a  source 


of  1  A-tum-m*.  Therefore,  if  a  source  magnetic 
irtomem  of  100  A-tum-m*  were  to  be  used,  then  the 
amplitude  of  the  magnetic  field  generated  should  be  2  x 
10^  A-cum/m.  Given  the  assumed  noise  level  of  the 
magnetic  field  of  2  x  10~*  A-tum/m.  the  resulting  ratio 
of  S/N  level  would  be  about  10.  Magnetic  moments  of 
SOO  A-tums  -m*  are  reasonably  easy  to  obuin,  thus 
developing  an  adequate  S/N  level  should  not  be  a 
problem. 

Figures  13  through  19  present  calculated  magnetic 
fields  at  the  se.isor  antenna  for  a  source  magnetic 
moment  of  100  A-tum-m*.  On  several  of  these  figures, 
the  magrthude  of  the  assumed  magnetic  noise  level,  2  x 
10~*  A-turn/m,  is  also  Indicated  for  comparison 
puiposes. 

Figiues  13a  and  b  indicate  that  differentiation 
between  sediment  types,  characterised  by  variations  in 
formation  factor  from  2.5  to  10.  is  possible  for  induction 
numbers  ranghtg  from  0.37  to  3  7.  This  range  in 
induction  numbers  corresponds  to  frequencies  from  42 
to  4200  Hs  (see  Equation  6) . 

System  10  m  above  sea  floor,  two-layer 
sea  floor,  deep  water 

Figures  14a  and  b  present  calculated  magnetic  fields 
at  the  seitsor  for  antennas  towed  10  m  above  the  sea 
floor  in  infinitely  deep  ocean  wdiere  the  sea  floor 
consists  of  a  6-m-thick  sediment  layer  over  an  infinitely 
thick  second  layer.  The  resistivity  of  the  top  layer  is 
input  as  0.62.  1.25.  and  2.50  dm  with  that  of  the 
bottom  layer  maintained  at  2.50  Hm  (Fig.  C-2  of 
Appendix  C).  As  found  for  the  previous  model,  the 
resistivity  variations  in  the  top  6-m-thick  layer  are 
clearly  observed  in  the  magnetic  field  measurements. 

It  is  significant  to  observe  in  Figures  14c  and  d  that 
the  magnetic  measurements  for  Case  1C.  a  low-resis¬ 
tivity  (toft?)  sediment  layer  of  infinite  thickness,  are 
difficult  to  distinguish  fiom  those  for  Case  2B, 
low-resistivity  (soft?)  layer  over  the  high-resistivity 
(stiff?)  sediment  layer.  These  dau  indicate  that  the 
surficial  6-m-thick  layer  has  a  dominating  influence 
over  the  magnetic  response  at  the  sensor,  and  that  the 
underlying  layer  2  sediments  will  be  more  difficult  to 
classify  by  electromagnetic  methods,  when  the  lower 
layer  has  a  somewhat  higher  electrical  resistivity. 

System  10  m  above  sea  floor,  water 
depth  20  m 

The  influence  of  shallow-water  depth  is  shown  in 
Figures  15a  and  b.  These  figures  present  results  for 
Cases  iA,  B  and  C,  Appendix  C,  which  are  the  same  as 
the  Cases  1  except  for  water  depth.  Again,  as  in  the 
previous  two  cases,  the  variations  in  the  sediment 
resistivity  from  0.6  to  1.25  to  2.50  Om  are  easily 
detected.  comparison  between  Figures  ISa  and  13a 
indicates  that  the  fields  measured  in  20  m  of  water,  as 
opposed  to  those  in  an  infinite  water  depth,  have 
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Figure  14(a).  Resolution  of  sediment  types,  array  10  m 
above  sea  floor,  6  m-thick  layer  1  over  thick  layer  2, 
real  component  of  hz  generated  by  HMD. 


Figure  14(b).  Resolution  of  sediment  types,  array  10  m 
above  sea  floor,  6-m  layer  1  over  thick  layer  2, 
imaginary  component  of  Hz  generated  by  HMD. 
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Figure  14(c).  Resolution  of  sediment  type  of  underly¬ 
ing  layers,  array  10  m  above  seafloor,  6-m  layer  1  over 
thick  layer  2,  real  component  of  Hz  generated  by 
HMD. 
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Figure  14(d).  Resolution  of  sediment  type  of  underly¬ 
ing  layers,  array  10  m  above  seafloor,  6-m  layer  over 
thick  layer  2,  imaginary  component  of  Hz  generated 
by  HMD. 
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Figure  15(a).  Resolution  of  sediment  types  in  20-m  I  Figure  15(b).  Resolution  of  sediment  types  in  20~m 

water  depth,  array  10  m  above  seafloor,  very  thick  j  water  depth,  array  10  m  above  seafloor,  very  thick 

•  layer,  real  component  of  Hz  generated  by  HMD.  layer,  imagirutry  component  of  Hz  generated  by  HMD. 
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Figure  16(a).  Irfluence  of  resistive  boundaries,  real  Figure  16(b).  lr\fluence  resistive  boundaries,  imaginary 
component  of  Hz  generated  by  HMD.  component  of  Hz  generated  by  HMD. 
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Figure  17(a).  Ir\fluence  of  distance  to  sea  surf  ace,  array  Figure  17(b).  Irtfluence  of  distance  to  sea  surface, 

10  m  above  sea  floor,  q=2.5  ohm-m,  real  component  array  10  m  above  seafloor,  c  =  2.5  ohm-m,  imaginary 

of  Hz  generated  by  HMD.  component  of  Hz  generated  by  HMD. 
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Figure  18(a).  Decrease  in  response  with  increasing  Figure  18(b),  Decrease  in  response  with  increasing 
distance  above  sea  floor,  q=2.5  ohm-m,  real  compo-  distance  above  sea  floor,  q^2.5  ohm-m,  imaginary 
nent  of  Hz  generated  ly  HMD.  component  of  Hz  generated  by  HMD. 
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INDUCTION  NUkIBER 

Figure  19(a).  Influence  of  uncertainty  in  height  of 
array  above  sea  floor,  q—2.5  ohm-m,  real  component 
of  Hz  generated  by  HMD. 

decreased  by  a  factor  of  two,  bjt  the  resolution  between 
resistivities  is  similar.  In  addition,  the  magnetic  response 
from  Cases  5  and  1  are  nearly  180**  out  of  phase,  with 
this  difference  caused  by  the  high  electrical  resistivity 
contrast  of  the  air>water  interface. 

Distortions  due  to  seafloor  and  ocean  surface 
Figures  16a  and  b  show  the  relative  magnitudes  of 
the  field  distortions  introduced  by  the  sea  floor  and 
ocean  surface  boundaries  considered  in  the  cases 
treated  above.  Case  lA  places  the  EM  system  in  an 
infinitely  deep  ocean  10  m  over  a  2.5  Hm  sedimentary 
half-space.  Case  5A  lowers  the  ocean  surface  to  just  10 
m  over  the  EM  system,  and  Case  6A  leaves  the  ocean 
surface  at  10  m  over  the  EM  system  but  lowers  the  sea 
floor  to  an  infinite  distance  away.  The  proximity  of  the 
ocean  surface,  a  resistive  boundary,  and  of  the  sediment 
surface,  also  a  resistive  boundary  as  compared  to  the 
water  column,  produces  reciprocal  phase  differences 
(Fig.  16a  and  b).  Note  that  within  the  accuracy  of  the 
simulation,  the  sum  of  the  magnetic  responses  from 
Cases  1  and  6  is  equal  to  the  response  calculated  for 
Case  S. 

Influence  of  the  distance  to  ocean  surface 
Figures  17a  and  b  indicate  the  relative  distortion  in 
signal  caused  by  towing  the  EM  system  through  an 
ocean  with  various  water  depths.  In  all  cases  considered, 
the  EM  system  is  shown  deployed  10  m  above  a  2.3  flm 
sedimentary  half-space.  The  depth  of  water  is  varied 
such  that  the  distance  to  the  ocean  surface  from  the  EM 
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Figure  19(b).  Influence  of  uncertainty  in  height  of 
array  above  sea  floor,  q=2.5  ohm-m,  imaginary  com¬ 
ponent  of  Hz  generated  by  HMD. 

system  is  infinity  (Case  1),  10  m  (Case  S),  and  5  m 
(Case  4).  The  results  displayed  in  Figures  17a  and  b 
show  that,  by  varying  the  towing  depth  of  the  EM  system 
belov:  the  ocean  surface  from  10  to  S  m,  a  variation  in 
magnetic  field  of  nearly  two  orders  of  magnitude  is 
induced.  This  large  influence  of  the  sea  surface  on  the 
measurement  indicates  that  the  EM  system  should  be 
deployed  at  considerable  depth/distance  from  the  sea 
surface  when  possible.  The  electromagnetic  field 
response  will  depend  heavily  on  any  resistive  bounda¬ 
ries.  Therefore,  it  follows  that  EM  measurements  should 
be  acquired  near  the  boundary  of  interest  (the 
seafloor/sediments)  and  away  from  a  high-contrast 
boundary  that  is  already  defined  (the  ocean  surface) . 


Influence  of  height  above  sea  floor 
Figures  18a  and  b  illustrate  the  variations  in  the 
magnetic  field  response  induced  by  varying  the  height  of 
the  EM  system  above  the  sea  floor.  In  all  cases,  the  EM 
system  is  in  an  ocean  of  infinite  depth  over  a  2.S  Om 
sedimentary  half-space.  The  distance  of  the  EM  system 
above  the  sea  floor  is  varied  from  10  m  (Case  lA)  to  20 
m  (Case  7A)  and  to  30  m  (Case  7fi).  By  raising  the  EM 
system  from  10  to  20  m  above  the  sea  floor,  the 
magnetic  response  measured  from  the  sediments 
decreases  by  an  order  of  magnitude,  with  an  equal 
change  in  raising  the  system  from  20  to  30  m.  This  large 
decrease  in  response  is  caused  by  the  large  electromag¬ 
netic  propagation  attenuation  induced  by  the  seawater 
column.  Note  that  the  magnetic  response  for  increasing 
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standoff  distance  is  not  only  attenuated,  but  the 
response  is  also  confined  to  a  lower  frequency  range 
(see  Eq.  4). 

Influence  of  elevation  errors 

Figures  19a  and  b  show  the  variations  in  the  vertical 
magnetic  field  at  the  sensor  induced  by  vertical 

variation  in  the  position  of  the  horizontal  dipole  EM 
system,  i.e.,  standoff  distances  of  9,  10,  and  11m  above 
the  sea  floor.  The  variations  in  magnetic  field  are  on  the 
order  of  35  x  ItT*  to  65  x  10'*  A-tums/m.  In 
comparison,  the  variations  in  magnetic  field  arising  due 
to  change  in  sediment  resistivity  from  0.60  to  2.50  flm 
for  a  fixed  standoff  distance  of  10  m  (Fig.  13a  and  b) 
vary  over  a  range  of  25  x  10‘®  to  35  x  10^ 
A-tums/m—the  same  order  of  magnitude  as  the  error 
caused  by  a  +1  m  variation/error  in  the  EM  system 
standoff  distance.  To  achieve  required  measurement 
accuracy,  the  standoff  distance  from  the  sea  floor  will 
have  to  be  known  to  better  than  0.1  m.  This  standoff 
distance  measurement  accuracy  may  be  provided  by  the 
combined  use  of  an  acoustic  system  in  conjunction  with 
the  interpretation  of  the  high-frequency  components  of 
thici  EM  data. 

Influence  of  transmitter/receiver  tilt. 

The  error  in  magnetic  field  measurement  caused  by 
tilt  (or  pitch)  of  the  transmitter/receiver  system  in 
airborne  systems  has  been  studied  in  detail  by  Son 
(1965).  His  work  indicated  that  error  due  to  tilt  could 
be  as  large  as  several  percent.  Preliminary  calculations 
done  in  this  study  for  a  submerged  system  indicate  that 
the  tilt  angle  must  be  measured  and  appropriate 
corrections  made  to  the  field  data  before  a  reliable 
interpretation  is  possible. 

Discussion 

The  preliminary  evaluation  of  EM  system  perform¬ 
ance  in  the  classification  of  sea-floor  sediments  shows 
that  the  concept  is  technically  feasible.  Implementation 
of  the  concept  will  not  be  without  problems.  Proper 
performance  of  the  system  requires  that  thie  trans¬ 
mitter/receiver  antennas  be  towed  within  10  m  of  the 
sea  floor,  that  the  actual  standoff  distance  above  the  sea 
floor  be  known  to  the  nearest  0. 1  m  and  possibly  slightly 
better,  and  that  the  pitch  of  the  antennas  be  measured 
to  the  nearest  O.!". 

Possible  design 

The  firm  design  of  a  system  would  be  premature  at 
this  point;  however,  a  short  summary  of  the 
requirements  for  an  underwater  electromagnetic  system 
is  provided.  A  frequency  domain  system  should  have  a 
rigid  support  connecting  the  transmitter  and  receiver  at 


a  flxed  separation  of  5  to  10  m.  The  source  and 
receivers  should  have  a  bandwidth  of  10-1000  Hz.  Due 
to  the  small  amplitude  of  the  expected  electric  fields 
(pV/m),  induction  coil  sensors  should  be  used  for 
magnetic  field  measurements.  The  transmitter  may  be 
short  electric  dipoles  with  lengths  of  1  to  2  m  or 
induction  coils  with  a  diameter  on  the  order  of  1  m.  The 
ship’s  power  supply  could  be  used  to  drive  the 
ttansmitter  by  converting  the  power  into  the  desired 
waveform  using  control  amplifiers.  A  shipboard  AC 
power  supply  of  208  VAC  is  capable  of  providing  a 
SO-A  current  source  to  either  an  insulated  loop  source 
or  a  short  electric  dipole.  This  would  provide  a  magnetic 
moment  of  500  to  1000  A-tums-m*  or  an  electric 
dipole  moment  of  50  to  100  A-m.  A  loop  with  50  to 
100  turns  of  No.  8  AWG  copper  conductor  with 
0.045-inch  polyethylene  insulation  would  be  adequate 
to  handle  these  current  levels  with  a  resulting  weight  of 
less  than  500  lb. 

Induced  polarization  (IP) 
methods 

In  both  the  DC  and  EM  methods  presented,  the 
measured  elecuical  conductivity  of  the  sea-floor 
materials  has  been  assumed  to  be  independent  of  the 
harmonic  content  of  the  transmitted  signal.  In  reality 
electrical  properties  do  depend  on  frequency  to  some 
extent,  but  the  degree  of  dependence  is  small  enough 
that  the  effect  can  be  ignored  when  applying  the  DC  and 
EM  methoQS  to  the  determination  of  sediment  layer 
electrical  resistivities  and  thicknesses.  On  the  other 
hand,  it  is  possible  to  design  measurement  methods  that 
utilize  this  frequency-dependent  property  (Xeller  and 
Frischknecht.  1966)  to  detect  the  existence  of  clay 
minerals  in  a  sediment  and  possibly  even  determine  the 
percentage  of  clay  minerals  in  the  sediment  (Wynn, 
1987).  The  method  is  not  new  (i.e.,  early  work 
described  by  Schiumberger;  1920);  however,  modem 
development  originates  largely  with  Bleil  (1953). 

The  electrical  polarization  most  conductors  acquire  in 
an  electric  field  disappears  instantaneously  when  the 
field  is  switched  off.  Some  materials,  however,  retain 
part  of  their  polarization  for  a  time  of  the  order  of  a  few 
hundred  milliseconds  and,  although  the  electrochemical 
mechanism  of  the  effect  is  not  fully  understood,  it  is 
associated  with  conductors  containing  disseminated 
metallic  ore  minerals  or  colloidal  particles  (clay 
minerals)  (Griffiths  and  King,  1965).  When  the  current 
flow  is  abruptly  interrupted,  the  electric  field  will  decay 
as  shown  in  Figure  20,  with  a  quick  drop  to  a  value  Vs, 
which  is  usually  a  small  fraction  of  the  static  field,  Vp. 
After  this  initial  drop,  the  field  will  decay  slowly  with 
measurable  size  for  periods  of  minutes  to  hours  (Keller 
and  Frischknecht,  1966).  The  depolarization  or  decay 
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curve  it  feneraUy  logarithmic  in  aiiipe,  although 
contributing  components  with  different  time  consunts 
can  usually  be  identified  (Sumner,  1979). 

Prior  to  1950,  all  IP  measurements  were  of  the 
time 'domain  type  that  used  the  waveforms  shown  in 
Figure  21.  A  simple  on-off  step-function  current  (Fig. 
21a)  was  inqjressed  in  the  earth  by  the  grounded 
contacts  and,  in  polarizable  materials,  a  modified 
waveform  similar  to  F.gure  21b  would  be  observed. 
Analysis  of  the  voluge  waveform,  i.e.,  the  ratio  of 
secondary  potential  (Vs)  to  maximum  primary  potential 
(Vp)  and  the  time-history  of  the  transient  voltcge  (V,), 
gives  a  measure  of  the  IP  effect  (Sumner,  1979). 


Figure  21 .  The  time-domain  transmitted  and  received 
voltage  waveforms,  showing  the  inducing  primary  cur¬ 
rent  I  being  detected  as  a  maximum  primary  voltage 
Vp.  when  current  is  turned  off,  voltage  drops  to  a 
secondary  level  V,  and  the  transient  voltage  V,  decays 
with  time.  (After  Sumner,  1979) 


After  1950,  frequency-domain  field  equipment  has  * 
played  an  increasing  role  in  the  measurement  of  IP  data. 
The  voluge  decay  (V()  inducted  by  the  IP  phenomenon 
is  time-dependent;  therefore,  the  resistivity  must  also 
have  a  frequency  dependence  since  the  time  domain 
dau  are  related  to  the  frequency  domain  measurements 
through  the  Fourier  transform.  Thus,  polarized 
materials  can  be  viewed  as  having  an  impedance  that  is 
frequency  dependent.  This  impedance  may  be 
measured  using  frequency  domain  techniques;  an 
example  of  the  current  and  voluge  waveform  patterns 
for  this  type  of  IP  measurement  is  shown  in  Figure  22. 
The  difference  in  apparent  resistivities  or  measured 
voluges  at  high  and  low  frequencies  is  used  to  calculate 
parameters  useful  in  estimating  the  degiii^e  of 
polarization  (Sumner,  1979). 

A  practical  way  of  implementing  this  approach 
requires  time  synchronization  of  the  transmitter  and  the 
receiver  to  permit  the  measurement  of  the  relati"~  phase 
difference  between  transmitted  and  receiv .  signals 
(Fig.  23).  The  phase  difference,  determineo  cither  in 
time  or  as  a  phase  angle,  gives  the  polarization  of  the 
intervening  earth.  The  resulting  harmonic  responses  can 
be  decomposed  into  in-phase  and  out-of-phase 
components  over  the  range  of  observed  frequencies, 
v^ch  may  in  turn  be  plotted  in  the  complex  plane  as  a 
series  of  routing  vectors.  Dau  in  this  form  are  known  as 
complex  resistivity  measurements  (Fig.  24)  (Sunmer, 
1979).  Polarizable  materials  exhibit  characteristic 
frequency  rtsponse/i,  with  the  frequency  response  of  the 
imaginary  component  being  diagnostic  of  i^eral  type. 
Further,  the  magnitude  of  the  maximum  imaginary 
component  U  diagnostic  of  the  percenuge  content  of 
the  mineral  type  (Wynn,  19(7).  Wynn  has  successfully 
used  the  IP  method  in  shallow  water  off  the  Georgia 
coast  to  evaluate  the  technical  feasibility  of  mapping 
ilemenite  (tiunium  ore)  deposits. 

In  addition  to  meullic  minerals,  fine-grained  c’ays 
and  fibrous  minerals,  such  as  serpentinite,  can  produce 
a  moderately  strong  IP  effect.  In  many  areas,  the  clay 
minerals  provide  a  polarization  response,  which  may  be 
easily  measured  with  conventional  IP  instruments. 

Induced  polarization  at  low  frequencies  cannot  be 
explained  in  terms  of  atomic  or  molecular  structure. 
The  imporunt  phenomena  contributing  to  IP  appear  to 
be  electrode  reactions  between  electrolytes  and  meullic 
mineral  grains,  variation  of  electrolyte  mobility  in 
fine-grained  materials,  and  related  phenomena  (Keller 
and  Frischknecht,  1966).  At  this  time,  there  are  several 
proposed  electrical  theories,  none  of  which  has  a  firm 
foundation  based  on  the  fundamenul  electrochemical 
nature  of  polarizable  materials.  Further,  electrochemis¬ 
try  holds  the  clues  to  the  reality  and  meaning  of  mineral 
discrimination  as  gained  from  spectral  IP  data.  Much 
more  research  remains  to  be  accomplished  before  the  IP 
mechanism  is  fully  understood  (Sumner,  1979). 
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Figure  22.  Frequency-method  waveforms,  showing  a  con¬ 
trolled  constant  inducing  current  I  at  frequencies 
and  being  detected  as  voltages  and  where 
V^<  V^.  The  dashed  line  is  the  sinusoidal  filtered 
voltage.  (After  Sumner,  1979) 

At  this  time  then,  it  is  established  that  some  clay 
minerals,  specifically  the  montmorillonites,  do  exhibit  a 
strong  IP  effect,  and  that  quartz  sands  exhibit  virtually 
no  IP  effect,  leading  to  the  hypothesis  that  we  may  be 
able  to  distinguish  cohesive  from  cohesionless  sediments 
by  using  IP  (Wynn,  1987).  The  fundamental 
electrochemical  nature  of  the  clay  minerals  is  not 
sufficiently  well  understood  to  support  the  development 
of  electrical  theories  for  the  prediction  of  the  IP 
re^nse  of  a  pure  chy  mineral  assemblage,  much  less 
the  usual  mixture  of  clay  minerals  and  other  sediment 
components.  Thus,  for  the  near>term,  a  clay  mineral 
classification  and  quantification  scheme  will  have  to  be 
empincal  in  nature  based  on  measured  data.  Further,  it 
will  probably  be  necessary  to  build  that  empirical  dau 
base  on  undisturbed  sediment  samples  and  possibly  on 
in  situ  measurements.  Slight  drying  of  geologic  materials 
is  known  to  alter  their  IP  response,  and  hydrostatic 
pressure  release  and  slight  sample  shear  deformation 
during  coring,  core  extrusion,  and  sample  preparation 
might  also  have  adverse  effea  on  the  IP  re^nse. 
These  potential  problems  can  all  be  circumvented,  and 
they  do  not  negate  the  idea  of  using  IP  for  sediment 
classification. 

IP  measurements  require  high-quality  antennas, 
transmitter,  and  receiver.  These  components  are  the 
same  as  those  used  in  a  high-quality  DC  resistivity 
system.  Thus,  the  same  hardware  and  electronics  can  be 
used  for  both  IP  and  DC  resistivity  (although  cheaper, 
lesser-quality  systems  can  suffice  for  DC  resistivity 
measurements). 

Overall,  then,  IP  is  a  very  attractive  technology  as  an 
add-on  to  DC,  EM,  and  acoustic  technologies  for 
near-surface  sediment  classification. 


Summaiy  and  conclusions 

Field  results 

•  Field  results  show  that  a  practical  DC  resistivity 
amy  of  40  m  length  dragged  on  the  sea  floor  should 
give  acceptable  resolution  of  subseafloor  layer  parame¬ 
ters  to  subseafloor  penetrations  of  6  m. 
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Analytical  results 

•  To  resolve  properties  of  a  thin  surficial  sediment 
layer  (1  m  thickness)  the  DC  array  must  be  deployed  on 
or  less  than  1  m  above  the  sea  floor. 

a  If  averaging  of  properties  over  a  thicker  layer  of  6  m 
thickness  is  acceptable,  then  resolution  can  be  achieved 
with  a  DC  array  standoff  distance  as  much  as  10  m. 

•  An  EM  array  is  predicted  to  resolve  resistivity 
variations  (and  sediment  type)  in  the  upper  6  m  of  the 
sea-floor  sediments  over  an  EM  bandwidth  of  several 
kilohertz  when  towed  within  10  m  of  the  sea  floor. 

•  The  EM  field  response  is  strongly  influenced  by  the 
air-sea  interface;  therefore,  the  EM  array  must  be 
operated/towed  far  away  from  the  sea  surface  and  close 
to  the  interface  of  interest,  the  sea  floor,  to  enhance 
resolution. 

e  Unknown  variation  or  error  in  the  height  of  the  EM 
array  above  the  sea  floor  causes  a  significant  error  in  the 
sediment  resistivities  calculated.  To  achieve  required 
measurement  accuracy,  the  distance  of  the  EM  array 
above  the  sea  floor  will  have  to  be  known  to  better  than 
0.1  m.  This  altitude  accuracy  is  achievable  through  the 
joint  use  of  acoustic  height  information  and  interpreta¬ 
tion  of  high-frequency  EM  data. 

e  Tilt  of  EM  system  coil  antennas  can  result  in  error 
of  predictions  of  several  percent:  thus,  tilt  angle  of  the 
antennas  must  be  measured  and  corrections  made  to  the 
data  to  achieve  a  reliable  interpreution.  Tilt  sensors  are 
available  to  provide  the  required  measurement  accuracy 
or  better. 

e  IP  technology  shows  promise  for  the  detection  of 
clay  minerab  in  surficial  sediments  and  possibly  even  for 
the  coarse  estimation  of  proportion  of  clay. 

Recommendations 

e  Field  data  should  be  obtained  to  validate  the 
computer  predictions  of  DC  resistivity  system  perform¬ 
ance.  A  DC  inverted  Schlumberger  array  capable  of 
being  towed  on  or  near  the  sea  floor  and  capable  of  6-m 
penetration  (about  40  m  long)  should  be  designed  and 
fabricated.  Deckside  components  should  be  borrowed 
or  rented.  Dau  should  be  gathered  over  area(s)  with 
adequate  ground  truth,  and  the  quality  of  the  data 
gathered  and  the  resolution  of  the  sediment  resistivities 
calculated  should  be  evaluated  and  compared  to  ground 
truth  data.  Electrode  noise  should  measured  in  these 
open  ocean  conditions  and,  if  necessary,  options 
developed  to  reduce  noise  levels  and  thus  improve 
signal-to-noise  levels  and  dau  resolution. 

•  The  initial  work  on  determining  EM  resistivity 
system  performance  should  be  extended  to  compare 
insulated  coil  transmitter  performance  with  that  of 
electric  dipole  sources,  i.e.,  additional  electric  and 
magnetic  dipole  resportses  should  be  calculated.  Initial 


work  indicates  that  electric  dipole  sources  may  provide 
adequate  fields  at  lower  cost  and  weight  than  insulated 
coil  transmitters. 

•  Time  domain  EM  techniques  should  be  evaluated 
for  application  to  the  surficial  sea-floor  resistivity 
measurement  problem  because  operating  in  the  time 
domain  may  simplify  hardware  design  requirements. 
Time  domain  systems  measure  the  secondary  fields  after 
the  primary  fields  are  terminated:  therefore,  time 
domain  systems  are  less  seitsitive  to  the  relative  position 
and  orienution  of  the  receiver.  Because  dimensionless 
tolerances  are  less,  the  transmitter  and  receiver  may  not 
require  a  rigid  boom  support. 

•  A  comparison  of  the  DC  resistivity  and 
electromagnetic  responses  over  a  suite  of  sedimentary 
models  should  be  done.  This  would  provide  a  measure 
of  the  relative  sensitivity  of  the  two  techniques  to 
variations  in  sea-floor  sediments.  The  ability  of  the  two 
methods  to  complement  each  other  and  perhaps  resolve 
certain  ambiguous  cases  should  be  considered  and 
evaluated. 

•  A  laboratory  and  fleld  IP  effort  should  be 
undertaken  to  determine  the  feasibility  of  detecting  clay 
mineral  existence  and  percent  content  to  a  degree  that 
would  provide  practical  input  to  remote  sediment 
classification. 

References 

Alfano,  L.  (1962).  Geoelectrical  Prospecting  with 
Underground  Electrodes.  Geophysical  Prospecting  10 
(3)  271-289. 

Archie,  G.  E.  (1942).  The  Electrical  Resistivity  Log 
as  an  Aid  in  Determining  Some  Reservoir  Characteris¬ 
tics.  Transactions,  AIME  146,  S4-62. 

Barnes,  B.  B.,  R.  F.  Corwin,  T.  G.  Hildenbrand,  L. 
Jackson,  R.  Kessler,  W.  Takeyama,  M.  Homick,  and  R. 
Jenkins  (1973).  Geophysics  Applied  to  Geotechnical 
Problems  in  a  Marine  Environment.  A  Case  Study: 
Monterey  Bay,  California.  U.S.  Department  of 
Commerce,  annual  report.  ARPA  Order  No.  1979, 
Amend.  3. 

Beard,  R.  M.  (1984).  Expendable  Doppler 
Penetrometer  for  Deep  Ocean  Sediment  Strength 
Measurements.  Naval  Civil  Engineering  Laboratory, 
Port  Hueneme,  California,  TR-905. 

Beck-.-in,  H.,  and  G.  Demiray  (1976).  The  Nanolog, 
a  New  Way  to  Log  the  Resistivity  of  the  Seafloor. 
Geophysical  Prospecting  22,  309-316. 

Bennett,  R.  H.,  D.  N.  Lambert,  M.  H.  Hulbert,  J.  T. 
Bums,  W.  B.  Sawyer  and  G.  L.  Freeland  (1983). 
Electrical  Resistivity/Conductivity  in  Seabed  Sediments. 
In  CRC  Handbook  of  Geophysical  Exploration  at  Sea, 
R.  A.  Geyer  and  J.  R.  Moore  (eds.),  CRC  Press,  Boca 
Raton,  Florida,  pp.  333-375. 


C* 


27 


Beyer,  J.  H.  (1V71).  Unpublished  repor:  reuined  by 
R.  F.  Corwin. 

Bischoff,  J.  (1978).  Die  Anwendung  Ceoelektrischtr 
Widersiandsmelhoden  in  Marinen  Bereich.  Ph.D.  thesis, 
Tech.  University  of  Berlin  (in  German). 

Bleil,  D.  F.  (1953).  Induced  Polarization:  A  Method 
of  Geophysical  Prosp>ecting.  Geophysics  18(3)  636. 

Bogolovski,  V.  A.  and  A.  A.  Ogilvy  (1974).  Detailed 
Electrometric  and  Thermometric  Observations  in 
Offshore  Areas.  Geophysical  Prospecting  22,  381-392. 

Boyce,  R.  E.  (1967).  Electrical  Resistivity  of  Modem 
Marine  Sediments  from  the  Bering  Sea.  M.  S.  thesis, 
San  Diego  Sute  College,  California. 

Chave,  A.  D.  (1983).  Numerical  Integration  of 
Related  Hankel  Transforms  by  Quadrature  and  Fraction 
Expansion.  Geophysics  48,  1671-1687. 

Coggon,  J.  H.  and  H.  F.  Morrison  (1970). 
Electromagnetic  Investigation  of  the  Sea  Floor. 
Geophysics  35,  467-489. 

Constable,  S.  C.,  C.  S.  Cox,  and  A.  D.  Chave 
(1986).  Offshore  Electromagnetic  Surveying  Tech¬ 
niques.  Expanded  Abstracts,  Fifty-Sixth  Annual  Inter¬ 
national  Meeting  of  the  Society  of  Exploration 
C-  ^physicists,  November  2-6,  Houston,  Texas,  pp. 
81-85. 

Corwin,  R.  F.  (1976).  Offshore  Use  of  the 
Self-Potential  Mnhod.  Geophysical  Prospecting  24(1), 
79-90. 

Corwin,  R.  F.  (1983).  Marine  Permafrost  Detection 
Using  Galvanic  Electrical  Resistivity  Methods.  Proceed¬ 
ings,  19r,  Offshore  Technology  Conference,  Houston, 
Texas,  Ui'C  4480,  pp.  329-336. 

Corwin,  R.  F.  (1985).  Electrical  Resistivity  Tech- 
niout  or  Offshore  Arctic  Geotechnical  Engineering 
App'i'.  ons.  Civil  Engineering  in  the  Arctic  Offshore. 
Proct  .igs,  Arctic  '85,  American  Society  of  Civil 
Engim  ,  pp.  137-143 

Corwin,  R.  F.  (1986).  Investigation  of  the  Direct- 
Current  Electrical  Resistivity  Method  for  Sea-Floor 
Property  C  ssification.  Report  prepared  by  Harding- 
Lawson  Assoc.,  Novato,  California,  for  the  Naval  Ocean 
Researc?'  •■d  Development  Activity,  NSTL,  Mississippi, 
under  Comract  No.  N00014-86-M-6035,  5  August,  p. 
39. 

Corwin,  R.  F.  and  U.  Conti  (1973).  A  Rugged 
Silver-Silver  Chloride  Electrode  for  Field  Use.  Review 
of  Scientific  Instruments  44,  708-711. 

Corwin,  R.  F.  and  U.  Conti  (1973a).  Sea-Floor 
Electrical  Resistivity  Measurement,  Using  a  Diver-Oper¬ 
ated  System.  Proceedings,  Offshore  Technology  Confer¬ 
ence,  1, 1-637  -  1-644. 

Francis,  T.  J.  G.  (1977).  Electrical  Prospecting  on 
the  Continental  Shelf.  Institute  of  Geological  Sciences. 
London,  Report  77/4. 


Fraser,  D.  C.  (1978).  Resistivity  Mapping  with  an 
Airborne  Multicoil  Electromagnetic  System.  Geophysics 
43.  144-172. 

Fraser,  D.  C.  (1979).  The  Multicoil  II  Airborne 
Electromagnetic  System.  Geophysics  44,  1367-1394. 

Fraser,  D.  C.  (1981).  Magnetic  Mapping  with  a 
Multicoil  Airborne  Electromagnetic  System.  Geophysics 
46.  1579-1593. 

Fullagar,  P.  K.  and  D.  W.  Oldenburg  (1984). 
Inversion  of  Horizontal  Loop  Electromagnetic  Fre¬ 
quency  Soundings.  Geophysics  49,  150-164. 

Giannini,  J.  and  D.  Thayer  (1982).  Extremely  Low 
Frequency  (ELF)  Quasi-Static  Propagation  Measure¬ 
ments  for  a  Calibrated  Electric  Field  Source  in  the 
Ocean.  L£££  Trartsactions  on  Antennas  and  Propaga¬ 
tion,  AP-30,  (5),  September,  pp.  825-831. 

Glenn,  W.  E.,  Jisoo  Ryu,  S.  H.  Ward,  W.  J.  Peoples 
and  R.  J.  Phillips  (1973).  The  Inversion  of  Vertical 
Magnetic  Dipole  Sounding  Data.  Geophysics  38(6) 
1109-1129. 

Griffiths,  D.  H.  and  R.  F.  King  (1965).  Applied 
Geophysics  for  Engineers  and  Geologists.  Pergamon 
Press,  New  York. 

Hutt,  J.  R..  and  J.  W.  Berg,  Jr.  (1968).  Thermal  and 
Electrical  Conductivities  of  Sandstone  Rocks  and  Ocean 
Sediments.  Geophysics  33(3),  489-500. 

Keller.  G.  V.  and  F.  C.  Frischknocht  (1966). 
Electrical  Methods  in  Geophysical  Prospecting.  Per¬ 
gamon  Press,  New  York. 

Kermabon,  A.,  C.  Gehin,  and  P.  Blavier  (1969).  A 
Deep-Sea  Electrical  Resistivity  Probe  for  Measuring 
Porosity  and  Density  of  Unconsolidated  Sediments. 
Geophysics  34(4),  554-571. 

Kong,  J.  A.  (1972).  Electromagnetic  Fields  due  to 
Dipole  Antennas  over  Stratified  Anisotropic  Media. 
Geophysics  37,  985-996. 

Lambert,  D.  N.  (1987).  An  Evaluation  of  a 

Computerized  Acoustic  Sediment  Classification  System. 
Naval  Ocean  Research  and  Development  Activity, 
NSTL,  Mississippi.  NORDA  Report  169  (in  press). 

Roy,  K.  K.  and  H.  M.  Elliot  (1981).  Some 

Observations  Regarding  Depth  of  Exploration  in  DC 
Electrical  Methods.  Geoexploration  19,  1-13. 

Schlumberger,  C.  (1920).  Etude  sur  la  Prospection 
Elect rique  du  Sous-Sol.  Gauthier- Villars  et  Cie.,  Paris. 

Schlumberger,  C.,  M.  Schlumberger,  and  E.  G. 

Leonardon  (1934).  Electrical  Exploration  of  Water- 
Covered  Areas.  Transactions,  American  Inst'tute  of 
Mining  and  Metallurgical  Engineers  110,  122-134. 

Son,  K.  H.  (1985).  Interpretation  of  Electromagnetic 
Dipole-Dipole  Frequency  Sounding  Data  over  a 
Vertically  Stratified  Earth.  Ph.D.  Thesis,  North 

Carolina  State  University. 


28 


Sumner.  J.  S.  (1979).  The  Induced-Polarization 
Exploration  Method.  In  Geophysics  and  Geochemistry 
in  the  Search  for  Metallic  Ores,  P.  J.  Hood  (ed.). 
Geological  Survey  of  Canada,  Economic  Geology  Rei>ort 
31,  pp.  123-133. 

Tang,  C.  M.  (1979).  Electromagnetic  Fields  due  to 
Dipole  Antennas  Em^dded  in  Stratified  Anisotropic 
Media.  IEEE  Transactions  on  Antennas  and  Propaga¬ 
tion  AP-27  (5).  665-670. 

Valent,  P.  J.  (1986).  Shaliow~Water  RecoverabU 
Penetrometer  System.  Proposal  for  assembly  to  Naval 
Oceanographic  Office  from  Naval  Ocean  Research  and 
Development  Activity.  NSTL,  Mississippi,  24  Novem¬ 
ber. 

Wait,  J.  R.  (1954).  Mutual  Coupling  of  Loops  Lying 
on  the  Ground.  Geophysics  19(2),  290-296. 


Weaver,  J.  T.  (1967).  The  Quasi-Static  Field  oi  an 
Electric  Dipole  Embedded  in  a  Two-Layer  Conducting 
Half-Space.  Canadian  Journal  of  Physics  45, 
1981-2002. 

Webb,  S.  C.,  S.  C.  Constable,  C.  S.  Cox,  and  T.  K. 
Deaton,  1985.  A  Seafloor  Electric  Field  Instrument,  J. 
Geomagnetics  and  Geoelectr.,  Vol.  37,  pp.  1115-1129. 

Won,  I.  J.  and  K.  Smits  (1985).  Airborne 
Electromagnetic  Bathymetry.  Naval  Ocean  Research 
and  Development  Activity,  NSTL,  Mississippi,  NORDA 
Report  94. 

Wynn,  J.  (1987).  Personal  Communication.  Office  of 
Mineral  Resources,  U.  S.  Geological  Survey,  Reston, 
Virginia,  8  January. 


29 


Appendix  A.  Calculated  resistivity  resolutions  for  DC 
Schlumberger  array 
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Figure  A-2.  Case  2,  DC  Resistivity  Resolution,  10m  Water  Depth, 
Array  on  Sea  Floor,  6m  Thick  Layer  1 
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Figure  A-3.  Case  3,  DC  Resistivity  Resolution,  lOm  Water  Depth, 
Array  on  Sea  Floor,  60m  Thick  Layer  1 
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Figure  A-4.  Case  4,  DC  Resistivity  Resolution,  10m  Water  Depth, 
Array  Im  Above  Sea  Floor,  im  Thick  Layer  1 
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Figure  A-7.  Case  7,  DC  Resistivity  Resolution,  Array  at  Sea  Surface 
in  10m  Water  Depth,  Im  Thick  Layer  1 
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Figure  A-8.  Case  8,  DC  Resistivity  Resolution,  Array  at  Sea  Surface 
in  10m  Water  Depth,  6m  Thick  Layer  1 
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Figure  A-12.  Case  12,  DC  Resistivity  Resolution,  Array  on  Sea  Floor 


In  100m  Water  Depth,  60m  Thick  Layer  1 
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Figure  A-19.  Case  19,  DC  Resistivity  Resolution,  Array  at  Sea  Surface 
In  100m  Water  Depth,  Im  Thick  Layer  1 
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Figure  A>20.  Case  20,  DC  Resistivity  Resolution,  Array  at  Sea  Surface 
In  100m  Water  Depth,  6m  Thick  Layer  1 
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Figure  A-21.  Case  21,  DC  Resistivity  Resolution,  Array  at  Sea  Surface 
In  100m  Water  Depth,  60m  Thick  Layer  1 
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Figure  A-22.  Case  22,  DC  Resistivity  Resolution,  Array  on  Sea  Floor 
In  1000m  Water  Depth,  Im  Thick  Layer  1 
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Figure  A-25.  Case  25,  DC  Resistivity  Resolution,  1000m  Water  Depth, 
Array  Im  Above  Sea  Floor,  Im  Thick  Layer  1 
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Figure  A-26.  Case  26,  DC  Resistivity  Resolution,  1000m  Water  Depth, 
Array  Im  Above  Sea  Floor,  6m  Thick  Layer  1 
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Figure  A-32.  Case  32,  DC  Resistivity  Resolution,  1000m  Water  Depth, 
Array  100m  Above  Sea  Floor,  6t,  Thick  Layer  1 
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Figure  A-33. 


Case  33,  DC  Resistivity  Resolution,  1000m  Water  Depth, 
Array  100m  Above  Sea  Floor,  60m  Thiel:  Layer  1 
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Appendix  B.  Analytical  model  for  transverse  electric  (TE) 
and  transverse  magnetic  (TM)  fields 
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The  electromagnetic  modeling  work  conducted  as  part  of  this  feasibility 
study  was  done  using  software  developed  at  NORDA.  This  modeling  package  was 
based  on  the  formulation  by  J.  A.  Kong  (1972)  and  generalized  for  buried 
dipoles  by  C-M  Tang  (1979).  The  algorithm  assumes  that  the  earth  is  com¬ 
posed  of  a  series  of  unia;.lal  anisotropic  layers.  The  details  of  the 
computational  scheme  are  the  same  as  those  presented  by  Tang  (1979)  with  the 
exception  of  the  derivation  of  tha  primary  field  terms.  There  exists  a 
discrepancy  in  the  calculation  of  these  terms  between  this  development  and 
the  earlier  derivation  by  Tang  (1979)  for  the  case  of  anisotropic  media. 

In  summary,  the  electromagnetic  fields  due  to  dipole  antennas  buried  in 
a  stratified  anisotropic  media  may  be  obtained  using  a  Transverse  Electric 
(TE)  and  Transverse  Magnetic  (TM)  mode  formulation.  The  technique  provides 
a  closed  form  solution  in  terms  of  a  set  of  Integral  equations  which  may  be 
numerically  evaluated  using  various  schemes,  Chave  (1983)* 

The  geometry  of  the  problem  is  indicated  in  Figure  B-1.  The  fields  are 
calculated  in  a  cylindrical  coordinate  system  (r,  0  ,  z)  with  the  z  axis 
directed  up.  A  time  dependence  of  exp(-lcot)  was  used.  The  sources  or 
various  electromagnetic  antennas  were  located  in  layer  (0)  with  m  layers 
above  and  n  layers  below  the  source  layer.  Each  of  these  layers  is  repre¬ 
sented  by  a  uniaxial  anisotropic  material,  with  the  conductivity  0  , 
permittivity  T  and  permeability  "JT  defined  by  the  following  tensors: 


0 

£' 

0 


0 

0 


Ci) 


The  electric  (E)  and  magnetic  (H)  fields  may  be  decomposed  into  two 
orthogonal  components,  TC  and  TM  modes.  The  field  components  from  each  of 
these  modes  may  have  vector  components,  one  of  which  propagates  parallel  to 
the  direction  of  the  layers  (transverse  components  (E.  ,  H,.))  and  the  other 

w  w 

normal  to  the  layer  interface  (vertical  component  (E  ,  H  )).  The  z  com- 

z  z 

ponents  of  the  fields  E  and  H  are  composed  of  only  TM  and  TE  modes, 
respectively.  These  components  may  be  defined  as  a  Sommerfeld  integral  as 
follows; 


(r)  =J  H,  dA 


When  the  z  components  of  the  fields  are  known,  the  transverse  components  may 
be  calculated  using  the  following  relations: 
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where 


V,-^  E^A) 

H,™(A)  =  (i!^)V,x(E^A)z) 
E,T^)  =  (-!^)  V,  X(H  JA)  i) 
H.^A)= 


The  governing  equations  for  the  electromagnetic  fields  In  a  region 
outside  of  any  source  are  provided  by  Maxwell’s  source  free  equations  which 
are  as  follows: 

VxT  =  icupH  Vx  H  =  -IcajIE 


These  equations  with  the  previously  defined  permeability  and  permit¬ 
tivity  tensors  may  be  transformea  into  two  sets  of  uncoupled  differential 
equations.  These  equations  may  then  be  solved  for  the  TE  and  TM  solutions. 

i 


In  the  case  of  the  TM  mode,  Maxwell's  equations  provide  a  differential 
equation  for  the  field  component.  This  has  the  following  form: 
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where 


which  may  be  transformed  as  follows. 


VfE,  +  =  0 


where 


^2  ^  -t-  yW  =  [A^a  -  ]  v. 

a  ^ 


The  solution  for  E  provides  the  TM  mode  response  for  the  general 

z 

electromagnetic  induction  problem. 


The  TE  mode  may  be  obtained  by  again  transforming  the  source  free 
Maxwell's  equations  into  a  differential  equation  which  this  time  is  a  func¬ 
tion  of  only  H  .  This  frequency  domain  relationship  has  the  following  form: 
z 


1  k2 

- ix - H 

b  az^  b 
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where 


-^1 

ar  ' 

[dri 

and  b 


JL- 


which  may  be  transformed  as  before  Into  the  following  form: 


V?H,  +  =  0 

where 

A»  =  l""*'  /"■>  =  [A^b  -  k2] 

b 

The  electromagnetic  respose  from  a  vertical  magnetic  dipole  or  small 
horizontal  loop  of  wire  with  a  time  varying  current  will  be  composed  of  only 
a  TE  mode  component,  when  the  source  is  embedded  in  a  layered  environment. 

An  example  o'"  this  situation  is  shown  in  Figure  B-2.  For  this  case,  a 
vertical  magnetic  dipole  is  submerged  in  the  ocean  with  a  sea  floor  composed 
of  horizontal  sedimentary  layers. 

The  response  generated  by  this  source-model  configuration  is  composed 
of  a  vertical  magnetic  field  a  radial  (horizontal)  magnetic  field  H,  and 
a  circular  (horizontal)  electric  field  Each  vector  field  is  composed  of 


a  primary  and  secondary  component. 


Hz  -  Hf  +  hJ 


H, E|  -eS+eJ 


Where  the  primary  fields  are  defined  as  follows: 


PRIMARY  FIELDS: 


,  AN  1(322  _  r2  |(2R2(R2_i2)  +  ^^^^2  _  3^2^ 
An  RS 

,  AN  rze**^  [(3  -  -  ik3R] 

471  R® 


\cou  IAN  re‘*^  m 
477  R3  ^ 
R  = 


-  IkRl 


2-h 


k  =  (w2^£  +  Icu^o]’'^*  £ 

<t)  =  2Tlf  ® 

M  =  PERMEABILITY  (Henry/m)  f 


a  PERMITTIVITY  (Farad /m) 
a  CONDUCTIVITY  (S/m) 
a  FREQUENCY  (Hz) 


The  primary  components  are  simply  those  field  components  which  would  be 
measured  if  the  source  layer  were  infinitely  thick.  The  secondary  field 
component  is  generated  by  variations  in  physical  parameters  at  each  layer 
interface. 
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The  total  field  response  for  an  arbitrary  number  of  layers  may  be 
formulated  as  a  set  of  Integral  equations.  The  TE  mode  response  is  given  by 
the  following  relations: 


°  ^  /  ^)  jj(A,)dA 

Jq 

“  A*  {  -  cj’  +  D™  }  J,(ir)dA 

“  f  +  cj’  W)#'?**  +  D™  }  J,Ur)cU 

Ja 


where 


ygn)  =  *  =  OBSERVATION  DEPTH 

(Ar),  J  »  0,1  BESSEL  FUNCTIONS  OF  ZERO  AND  FIRST  ORDER 

(A)  AND  D?^  (A)  COEFFICIENTS  DUE  TO  PHYSICAL  PROPERTY 

VARIATIONS 


The  first  term  under  the  integral  is  the  primary  field  component  and  the 
second  and  third  terms  represent  the  secondary  fields.  These  secondary 


components  may  easily  be  calculated  by  generating  the  and  coefficients 

and  then  numerically  integrating  the  resulting  kernel.  The  coefficients  are 
defined  as  follows: 


cj'  (A) 


Ant)  •ZrTdc 

Cl  e 


D«  (A) + , 


(A)  =  e 


2/0  d-i 


e 


■2/0  do 


(mi  (m). 

/oh  '/oh 

e  °  +  e  ° 


(m)  (m)  2y^ 
4.1  • 


(d,  -  do) 


where 


I  .  irra.  .te\ 
I  M,.,  -t-  «  '  ^  Q  ) 

( .  •  2rr<l« 

( 


+  e’ r’') 


.  .  •2r?’it,  pis) 

\1  +  M.1,0  ®  "  / 


d.,  -  DEPT/1  TO  BASE  OF  SOURCE  LAYER 
dg  B  DEPTH  TO  TOP  OF  SOURCE  LAYER 
h  B  DEPTH  TO  SOURCE 


These  coefficients  are  related  to  the  physical  properties  of  not  only  the 
source  layer  but  to  the  relative  property  variations  between  the  other 
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layers  as  well.  This  relationship  to  the  other  layers  is  indicated  ex- 


TE  TE 

plicitly  through  the  Q  and  R  terms.  These  terms  may  be  calculated  via 

TE 

the  following  recursive  relations  for  R  : 


i  ss  HI,  -n+l,  -n+2  •••  *2 


and  for  Q 


TE, 


qTE 


El 


n.l-1  \  1  + 

■’j+i 

0  j  =  m-1,  m-2,  m-3,  ...  1 


TE  TE 

When  the  R  and  Q  terms  have  been  generated  the  calculation  of  the  kernal 
function  is  a  straightforward  procedure.  The  remaining  integral  equation  is 


then  evaluated  using  an  adaptive  quadrature  scheme  in  conjunction  with  a 
continued  fraction  algorithm  which  easily  evaluates  the  slowly  convergent 
integrals,  Chave  (1983) • 

The  above  equations  may  be  transformed  into  the  TM  mode  solution  for 
the  E^,  field  components.  This  is  easily  done  by  making  the  follow¬ 

ing  substitutions: 


ej  =  e,  + 

w 


Vfi 


The  resulting  TM  field  response  is  generated  by  a  buried  vertical  electric 
dipole  in  a  layered  media. 


The  last  ease  to  be  considered  would  be  related  to  the  response  gener¬ 
ated  by  horizontal  electric  or  magnetic  dipoles.  The  resulting  secondary 
fields  generated  in  a  layered  environment  by  these  sources  have  the  follow¬ 
ing  components: 

H,,  E^,  Ep  E^ 
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The  evaluation  of  these  fields  proceeds  in  a  manner  similar  to  the 
previous  case  (TE  mode)  except  that  these  fields  have  a  much  more  complex 
kernel  function  which  has  a  contribution  from  both  TE  and  TM  modes. 


LAYER  m 


LAYER  m-1 
LAYER  m-2 


2 


LAYER  J _ 

LAYER  0 _ ^ 

LAYER  -1 _ 

LAYER  -2 


«>i 

".1 

'*2 


LAYER 
LAYER  -n 


FIGURE  B-1  rODEL  GEa-ETRY 
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°*>  OCEAN 


SEDIMENTS 

O  TRANSMIHER  (VERTICAL  MAGNETIC  DIPOLE) 

®  RECEIVER  {H„  H„  E^) 


•  FIGUf€  B-2  VERTICAL  mGNETIC  DIPOLE  SIBOGED  IN  THE  OCEAN 
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t  mA  lULAA J 


Appendix  C.  Sketches  of  seven  cases  examined  in  EM 
analytical  studies 


THETA 

(CIRCULAR) 


RADIAL 


WATER 


Py^  =  0.25  ohm-m 


r  =  10.  m 


TRANSMIHER 


-~w  »  1 

^  !  K  _  . 


RECEIVER 


h  =  10.  m 


SEDIMENT 


P^  =  2.5.1  1.25.»  AND  0.625  ohm-i 


Figure  C-1.  Cases  lA,  B,  and  C;  Resistivity  Study 


RADIAL 


r  =  10.  m 


Py^  =  0,25  ohm-m 

TRANSMIHER 


THETA 

(CIRCULAR) 


mJEH 


RECEIVER 


I  h  =  10.  m 


id  —  6.  m  P^-  1.25  and  0,625  ohm-m 


SEDIMENT  LAYER 


P2  =  2,5  ohm-ni 


SEDIMENT  BASEMENT 


Figure  C-2.  Cases  2A  and  B,  Resistivity  Study 
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^ - -  RADIAL 

THETA 

(CIRCULAR) 


ohm-m 

TRANSMITTER* 


T\ 


RECEIVER 
h  =  11.,  9.  m 


SEDIMENT 
Pi=  2.5 


ohm-m 


Figure  C-3.  Cases  3A  and  B,  Resistivity  Study 


Z 


- -  RADIAL 

THETA 

(CIRCULAR) 

AIR' 

^air  ^ 


r  =  10.  m 


d  =  15.  m 

TRANSMIHER 
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Figure  C-4.  Case  4A,  Resistivity  Study 
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AIR 

Pair  =  ** 


I - ^  RADIAL 

THETA 

(CIRCULAR) 


r  =  10.  m 


d  =  20.  m 

TRANSMIHER 
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i 

I 

I  h  = 

I 

i _ 


—RECEIVER 
10.  m 


WATER 

Pw  “  0,25  ohm-m 


SEDIMENT 

Pi  =  2.5.  1,25.  and  0.625  ohm-m 


Figure  C-5.  Cases  5A,  B,  and  C,  Resistivity  Study 
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Figure  C-6. 


Case  6A,  Resistivity  Study 
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THETA 
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r  =  10.  m 
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^  ©—RECEIVER 
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Py/  “  0.25  ohm-m 


h  =  20.,  30.  m 
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Pj  =  2,5  ohm-m 


RADIAL 


Figure  C-7.  Cases  7A  and  B 


Resistivity  Study 


UNCLASSIFIED _ 

SECURITY  ClASSinCATION  OP  THIS  PAGE 


11.  REPORT  SECURITY  CLASSinCATION 

Unclassified 


a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLA8SIFICAT10NIOOWNORAOINQ  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NORDA  Report  21 1 


6.  NAME  OF  PERFORMING  ORGANIZATION 


AD-A187  447 


lb.  RESTRICTIVE  MARKINGS 

None 


3.  DtSmiBUnON/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution  is 
unlimited. 


B.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

NORDA  Report  21 1 


71.  NAME  OF  MONITORING  ORGANIZATION 


Navai  Ocean  Research  and  Deveiopment  Activity  Navai  Ocean  Research  and  Development  Activity 


7b.  ADDRESS  fCIty,  SUM,  intf  ZIP  Cotfll 

Ocean  Science  Directorate 
NSTL,  Mississippi  39529-5004 


Sc.  ADDRESS  (aty,  SIM,  Md  ZIP  CodP) 

Ocean  Science  Directorate 
NSTL,  Mississippi  39529-5004 


U.  NAME  OF  FUNDING/SPONSORINQ  ORGANIZATION  Sb.  OFFICE  SYMBOL 

Naval  Ocean  Research  and  ottpmictbiti 

Development  Activity 


Sc.  ADDRESS  (City,  SIM,  tnd  ZIP  Cod*) 

Ocean  Science  Directorate 
NSTL,  Mississippi  39529-5004 


B.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 

PROJECT 

TASK 

ELEMENT  NO. 

NO. 

NO. 

62435N 

03585 

MOG 

II.  TITLE  fineludi  Stcurlty  Clttllllcttlon) 

Rapid  Underway  Sediment  Ciassification  by  Electrical  Methods 


12.  PERSONAL  AUTHOR(S) 

Philip  J.  Valent,  Edward  C.  Mozley,  and  Robert  P.  Corwin* 


ISl.  TYPE  OF  REPORT  13b.  TIME  COVERED  14.  DATE  OF  REPORT  (Yr.,  Mb.,  D»y) 

Final  I  From _ To _  _  August  1987 


IS.  SUPPLEMENTARY  NOTATION 

*  Consulting  Geotechnical  Engineer,  El  Cerrito,  California _ 


17.  COSATl  CODES  |  18.  SUBJECT  TERMS  (Contin't0  on  rovorto  ff  noeuury  fdontify  by  block  numbor) 

sediments,  electrical  conductivity,  direct  current  resistivity, 
electromagnetic,  induced  polarization 


IS.  PAGE  COUNT 

76 


19.  ABSTRACT  (Conllnu*  on  r*nrs*  II  ntetssiry  ind  Idtntify  by  block  numbtr) 

The  initial  results  of  a  technical  feasibility  study  for  using  direct  current  (DC)  electric  resistivity,  elec¬ 
tromagnetic  (EM),  and  induced  polarization  (P)  methods  to  classify  surficial  sea-floor  sediments  are 
reported  herein.  Analytical  work  performed  in  this  study  and  a  review  of  related  work  show  that  both 
the  DC  and  EM  resistivity  techniques  can  provide  electrical  measurements  of  sufficient  resolution  to 
classify  sediments  according  to  type  (mud,  clay,  silt,  sand,  gravel,  etc.).  Both  DC  and  EM  antennas 
must  be  close  to  the  sea-floor  surface  (within  10  m)  so  that  surficial  sediment  resistivities  can  be  suffi¬ 
ciently  resolved  to  permit  sediment  classification.  To  resolve  the  resistivity  of  a  1-m-thick  surficial  sea¬ 
floor  sediment  layer,  the  DC  antennas  must  be  deployed  less  than  1  m  above  the  sea  floor.  Both  DC 
and  EM  methods  can  adequately  resolve  the  resistivity  of  a  6-m-thick  surficial  layer  at  a  standoff  distance 
of  10  m.  Noise  tests  with  an  inverted  Schlumberger  array  towed  in  the  water  column  and  on  the  sea 
floor  show  noise  levels  to  be  within  acceptable  limits.  The  level  of  development  of  the  IP  method  for 
clay  mineral  detection  is  reported;  this  method  shows  considerable  promise  as  a  contributor  to  a  remote 
underway  sediment  classification  system. 


20.  DISTF.  (LABILITY  OF  ABSTRACT 

UNCLASSIFIEDiuhclMITED  □  SAME  AS  RPT 


22i  NAME  OF  RESPONSIBLE  INDIVIDUAL 


Philip  J.  Valent 


DD  FORM  1473,  83  APR 


21.  ABSTRACT  SECURITY  CUSSIFICATION 

DTic  USERS  □  Unclassified 


22b.  TELEPHONE  NUMBER  Unclad*  Ar*t  Cod*) 

(601)  688-4906 


EDITION  OF  I  JAN  73  IS  OBSOLETE. 


22c.  OFFICE  SYMBOL 

Code  363 


_ UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


